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PREFACE. 

f I^HE following lectures were collected from manuscripts left 
-*- by the late Professor J. C. Adams, and are now reprinted 
without change from his Collected Scientific Papers, Vol. II., 
pp. 1—84. 

It was thought that the wide interest attaching to the 
lunar problem reached many besides the professed astronomer, 
and would justify a separate publication of this short work. 

It is known that Adams contemplated the publication of 
some such essay himself, and it must be a matter of regret to 
all that he never did so. No pains have been spared to 
present the material properly, but it is unavoidable that it 
should appear from the hands of an editor in a less perfect 
form than if the author had issued it himself. 

Yet, allowing for this disadvantage, I think those best 
qualified to judge will consider this work fully worthy of 
Adams's great name. Of current elementary theories it may 
be said that they leave off where the difficulties of the subject 
begin, that is to say, where the various cases of slow con- 
vergence have been exposed, but not dealt with. It is perhaps 
not too much to say that these lectures carry us to the point 
where such difficulties end, in an adequate evaluation of all 
the chief constants. They leave the problem effectively solved 
and not merely stated, and shew the path clear for the for- 
mation of a detailed theory, if that is desired. 

R. A. SAMPSON. 

Durham. 
8 October, 1900. 
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LECTUBES ON THE LUNAR THEORY. 

[Lectures on the Lunar Theory were given by Adams 
from 1860 with few intermissions until 1889. Originally their 
aim was to illustrate geometrically the analytical processes and 
thereby render them more comprehensible, and they included 
some elegant theorems on the geometry of conies which have 
since become common property ; but every year several lectures 
were rewritten, and thus the whole fabric gradually changed 
into the form in which it is here presented, — the form, prac- 
tically, in which he gave them last. 

Perhaps it is superfluous to say that these Lectures stand 
upon a different footing to treatises that are intended to form 
the basis of Tables. With such, completeness is the first object 
and manner of presentation is secondary. Immense as is the 
labour of forming a treatise of this description, there exist 
several that leave little to desire in respect to fulness of 
detail. Indeed it may be suspected that their very perfection 
in the quality they profess has stifled to some degree the 
proper development of the subject, because at first sight it 
suggests that there is little left to do in the Lunar Theory, 
unless one is prepared to track down the inconsiderable errors 
that have eluded his Masters. This seems a mistake; the 
methods most suitable for the whole task adapt themselves 
comparatively ill to each detail of it, and there seems much 
that remains to be done in respect to inventing methods 
suitable for attacking separately, as far as they permit of 
separate attack, the many difficulties into which the theory 
divides at the outset, and frhence perhaps approximating to 
a. l. 1 



2 LECTURES ON THE LUNAR THEORY. 

a more adequate knowledge than we now possess of the relative 
motion of Three Bodies. So far, with the notable exception 
of Dr G. W. Hill and those that have followed him, we have 
seen comparatively little effort in this direction. 

This was the cardinal feature of Adams's plan, and his 
lectures shew the methods he had gradually elaborated to 
accomplish it. They separate the inequalities from one another 
as far as possible, and are content with indicating the manner 
in which these separate inequalities afterwards combine. To 
shew that, with so slight an apparatus and within so small 
a compass, the result is no mere sketch, we need but set 
side by side the coefficients of longitude found in these 
Lectures and the corresponding terms in Delaunay's TMorie. 

Adams. Delaunay. 



Variation, coeff. of 


sin 2D 




2106-4 


2106-25 




sin4D 




8-74 


8-75 


Parallactic inequality 


, sinD 




- 124-90* 


- 127-62 




sin 3D 




0-73 


0-84 




sin 5D 




001 


001 


Annual equation, 


sinZ' 




- 658-9 


- 659-23 




sin (2D - 


-V) 


152-09 


15211 




sin (2D + 1') 


- 21-57 


- 21-63 


Evection, 


sin (2D - 


-I) 


4596-6 


4607-77 




sin (2D + 


1751 


174-87 


Further, 










Motion of Apse, 


1- 


-c 


•0085, 


54 -008572 


Motion of Node, 


9- 


-1 


•003997 -003999 



For those to whom the difficulties of the Lunar Theory 
are known, these numbers need no comment. 

No Manuscript exists of Lecture L It is taken substantially 
from my own notes of 1889.] 

* With Delaunay's value of the Sun's Parallax, viz. 8"-75. 



LECTURE I. 



HISTOBICAL SKETCH. 



[The Lunar Theory may be said to have had its commence- 
ment with Newton. Many irregularities in the Moon's motion 
were known before his time, but it was he that first explained 
the cause of those irregularities and calculated their amounts 
from theory. 

Of the inequalities which are due to the action of the 
Sun, the first, — which is called the Evection, — was discovered 
by Ptolemy, who lived at Alexandria in the first half of the 
second century of our era, under the reigns of Hadrian and 
Antoninus Pius. At a very early period the relative distance 
of the Moon at different times could be told from the angle 
it subtended, and its orbit could thus be mapped out. By 
such means Ptolemy found that its form was not the same 
from month to month, and that the longer axis moved con- 
tinually though not uniformly in one direction. He represented 
this change by a motion of the centre of the ellipse, as we 
would put it, in an epicycle round the focus, obtaining thus 
a variable motion for the longer axis and a variable eccentricity. 

The . representation of position by means of epicycles is 
intimately related to the modern method of developing the 
coordinates in harmonic series; thus if we have 

00 = A 1 cos (n^ + a 1 ) + A 2 cos {nj, + ar 2 ) + • ■ ■ 

y = A 1 sin (njk + <*i) + A 2 sin (nj; + a 2 ) + . . . 

the motion of the point (%, y) is that on a circle of radius A 1 
with angular velocity n u around a centre which moves on a 

1—2 
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circle of radius A. 2 with angular velocity «.,, and so on; and 
if, more generally, we have 

x = A 1 cos (nj, + eti) + . . . 

y = B x sin (nj, + a^ + . . . 

we may reduce this case to the former by rewriting 

x = H (A\ + -Bi) cos (nj + a^ + „ (A^ - BJ cos (- nj-a^ + ... , 

y = 2 (A + Bl ^ sin ^ ht + a ^ + 2 ( Al ~ Sl ^ sin (~ 1ht ~ a i) + ■ • • • 

Probably we have here the reason why circular motions 
and epicycles were first employed. 

Tycho Brahe (1546 — 1601) discovered the existence of 
another inequality in the Moon's Longitude quite different 
from the Elliptic Inequality and the Evection. He found it 
bore reference to the position of the Sun with regard to the 
Moon ; so that when the Sun and the Moon were in conjunction 
or opposition or quadratures the position of the Moon was 
quite well represented by the existing theory, but from con- 
junction to the quadrature following, her position was more 
advanced than the place assigned to it, reaching a maximum 
of some 35' about half-way ; and in the second quadrant it 
was just as much behind. This inequality he called the 
Variation; it was the first that Newton accounted for theo- 
retically, and if we were to suppose the Moon and Sun to 
move, except for mutual disturbance, in pure circles in the 
same plane, it is the only one that would present itself. 

The next significant step was made by Horrox (1619 — 1641) 
who represented the Evection geometrically by motion in a 
variable ellipse, and gave very approximately the law of varia- 
tion of the eccentricity and the motion of the apse. He 
supposed the focus of the orbit to move in an epicycle about 
its mean place. 

Newton's Principia did not profess to be and was not 
intended for a complete exposition of the Lunar Theory. It 
was fragmentary; its object was to shew that the more 
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prominent irregularities admitted of explanation on his newly 
discovered theory of universal gravitation. He explained the 
Variation completely, and traced its effects in Radius Vector 
as well as in Longitude; and he also saw clearly that the 
change of eccentricity and motion of the apse that constitute 
the Evection could be explained on his principles, but he did 
not give the investigation in the Principia, even to the extent 
to which he had actually carried it. The approximations are 
more difficult in this case than in that of the Variation, and 
require to be carried further in order to furnish results of 
the same accuracy as had already been obtained by Horrox 
from observation. He was more successful in dealing with 
the motion of the node and the law of change of inclination. 
He shewed that when Sun and Node were in conjunction, 
then for nearly a month the Moon moved in a plane very 
approximately, and that the inclination of the orbit then 
reached its maximum, namely, 5° 17' about; but as the Sun 
moved away from the Node the latter also began to move, 
attaining its greatest rate when the separation was a quadrant, 
and that at this instant the inclination was 5° very nearly. 
He also assigned the law for intermediate positions. The 
fact that there was no motion when the Sun was at the 
Node, that is, in the plane of the Moon's orbit, confirmed 
his theory that these inequalities were due to the Sun's action. 

When we spoke of Newton's results as fragmentary and 
incomplete, let it not be understood that he gave only very 
rude approximations to the truth. His results are far more 
accurate than those arrived at in elementary works of the 
present day upon the subject. 

After Newton, Clairaut (1713 — 1765) treated the Lunar 
Theory analytically. He readily found the Variation and many 
other inequalities, but met with a difficulty in determining the 
motion of the apse. At first he made its mean motion only 
about one-half of the observed value, and supposed that this 
indicated a failure of Newton's law of the inverse square of 
the distance; but soon he recognized an error, caused by 
omission of terms which he had imagined would not affect 
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the result. When these were included the calculated amount 
was nearly doubled. 

The first Tables of the Moon which were sufficiently accu- 
rate for use in determining longitudes at sea by observation 
of Lunar Distances were those of Mayer. They obtained a 
prize offered by our Board of Longitude, and were published 
in 1770 by Maskelyne, the Astronomer Royal. 

The first Theories which gave the Moon's place with an 
accuracy equal to that of observation were those of Damoiseau 
and Plana. The former was published in 1827, preceded in 
1824 by Tables; the latter was published in 1832. 

Hansen's Tables, which are those now used, were constructed 
from theory and were published in 1857 at the expense of the 
British Government.] 



LECTURE II. 

ACCELERATIONS OF THE MOON RELATIVE TO THE 
EARTH. 

When three bodies move under their mutual attraction, 
their motions are unknown to us except in the cases when 
they are approximately elliptical ; but this restriction includes 
almost all the most important cases in the Solar System. 

If one body of the system is greatly predominant and if 
the lesser bodies are not close together, the centre of gravity 
of the greater body may be taken as a common focus around 
which the others move in approximate ellipses. Or again, if 
two bodies lie close together, their relative motion may be 
approximately the same as though they were isolated, although 
the system contains a third greatly predominant body; for 
their relative motion is affected by the difference of the attrac- 
tions of the central body upon them and not by the absolute 
value of those attractions. 

The Sun and Planets are an example of the first kind; 
the Earth, Moon and Sun of the second. The Earth and 
Moon describe orbits round the Sun which are approximately 
ellipses, and the Moon might be regarded as one of the planets; 
but this point of view would not be a simple one ; the disturb- 
ing action of the Earth would be too great, though it is never 
so great as the direct attraction of the Sun, that is to say, 
never great enough to make the Moon's path convex to the 
Sun. The more convenient method is to refer the motion of 
the Moon to the Earth, and counting only the difference of 
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the attractions of the Sun upon the Earth and upon the Moon, 
to find how this distorts the otherwise elliptical relative orbit. 
This is the method of the Lunar Theory. 

The position of the Sun must be referred to the same 
origin ; but since the Earth describes an ellipse about the Sun 
which is disturbed by the action of the Moon, if we choose as 
origin the Earth's centre, we must allow for the disturbance of 
the Sun's position by the Moon. This correction may be evaded 
by choosing as origin, not the Earth's centre, but the centre of 
gravity of the Earth and Moon, 
with respect to which the Sun 
describes a curve so closely el- 
liptical that no allowance is 
required. For, if S, E, M denote 
respectively the Sun, Earth, 
and Moon, and G the centre 
of gravity of E and M, the 
accelerating forces of S are 

on E S/SJS* in ES, 

oaf S/SM* in MS; 
and these imply accelerations of G of amount 

e^tm m* parallel t0 ES> 

ETMSM* P arallel t0 MS ' 
now the accelerations of S are 

E/SE* in BE, 
M/SM* in SM; 
hence the acceleration of G relative to S is 
S + E + M E .. . t „„ 

E + M S~W parallel t0 E8 > 
S+E+M M 
"E + M" SM' P arallel t0 MS > 



or 



/„ GE ., GM\ . _,„, 

e+m [ E -m- M 'm*) m0U > 

S + E + M f r SG , „ SG\ . __ 

-e+m- { E -sE> +M -m>) mGS - 
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Let EM = r, SG = r', SGM=m; then 

E + M ' 
Hence 



GE=-J^- Ji r. 

E + M 



E + Mr 



-, 3 cos to 



3 15 



M s r' 3 L 

/ E r\*(_ 3 
+ \E + Mrj { 2 + 

1 - 1 Ti M r 3 
E* r' 3 l E + Mr 



COS 2 ft) + 



)♦ ]■ 



iCOSO) 



/ H rV 



U' + ifr', 
and the accelerations of G are 



3 15 , 

■g + yCOS^ft) + 



T EM r 2 



COS ft) 



J 



in GiV, 



S + E 



7 + Jff", , £Af i*/ 3 15 ,\ 1 . _ 



Now r/r' is approximately ^x ; neglecting the square of this 
quantity, we see that S moves about G in a pure ellipse. 

Consider now the accelerations of the Moon relative to the 
Earth ; subtracting the accelerations of the Earth from those of 
the Moon, we find 

E + M 



ME' 



r a /MQ EG\. .._ 

+s [sm + sEi) mMG - 
s {sw* ~ m) parallel t0 GS; 



let E+M =fL, S = m ; then these become 
1 

E-Mr 



ft mr 



E-M r „ 
+ ^ — s> -/ o cos to 

E + Mr' 



mr 



3 cos a + -pj — ^ -f { — s + — cos 2 6) 
E + Mr'x 2 



3^15 

U 2 



1 in ME, 

)♦ ] 

parallel to GS. 
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In the accompanying spherical 
triangle, let G be the centre of the 
sphere, SM' the" ecliptic, and M' 
the projection of M. 

Let 1/u be the projection of 
ME on the plane of the ecliptic; 

6 the longitude of the Moon as seen from the Earth, 

6' the longitude of the Sun as seen from G, 

s the tangent of the Moon's latitude MM'. 

Then 

SM=a, 8M' = 6-ff, 

r = (1 + s^f- iir 1 , cos &> = cos (6 - 6') (1 + s 2 )-*, 

and the accelerations of M relative to E are 



jnt?_ m'(l+s 2 )* 
1 + s" r' 3 u 



m 



1 + 



E-M 1 
E + M r'u 



3 cos (6-6 r ) + 



in ME, 



r'u 



3 cos (6 - 6') + 



E-M 1 

E + M r'u 



.|(l +s 2 ) + ^cos 2 (0-0')) + . 



parallel to GS. 



Call these quantities U and V respectively; then if we 
resolve parallel to M'G, perpendicular to M'G in the plane of 
the ecliptic, and perpendicular to the plane of the ecliptic, we 
have the following quantities which we call P, T, 8 ; viz. : — 

P= [7(l+s 2 )-*-Fcos(6>-0'). 

T=-Vsm(6-6'), 

S= Us(l+s 2 )~i; 
and also 



S-Ps = 


Vs cos (6-6'). 


From these we find 


/am 2 m! 
(l+sf r' 3 u 


1 3 n/a ^ E-M 1 

S + jr COS 2(0 - #') + ,=, ry — 

L2 2 v ; E + M r'u 


{(S-I 


A cos (0 - 0') + y cos 3(61 - 6>')1 + ... 
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T=- 



r'u 



S-Ps= ^ 
r s u, 



s sin 2(0- 6') + -^ — .-, — 
2 v ' E + Mr'u 

- | s^j sin {6 -0')+ 1 ^ sin 3 (0- 0')| + . 



"3 3 _,„ „E-M 1 



{(f-l s2 ) cos ^-^ )+ ¥ cos3 ^-^} + ' 

Hence with the time as independent variable we have the 
equations of motion 

dP \dt) ~ r > 

rdt\ dt) 

d* 



dt:- 



(rs) = -S. 



Or we may write these with 6 as independent variable ; let 

dd 
dt' 

dd 

dt 



so that 
Then 






whence 



rrdH = T 

dd ~u 3 ' 

d¥ _ h-2 a &U ^du u dH 
dP t±U W U dd n dd' 



again, 



dt 
whence 



, s tt„ J d*s d^iA , rjdH .( ds du\ 

( rs ^ H ' u \ u d^- s d^ +H de u [ u d0- s dd)' 



H 



■«•(: 



d*s 



dH ds 



d 6* + s ) + H d8 ul de = Fs - s > 
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or the equations of motion may be written 



de ~u 3 ' 



Hv ^ + S )- Ps - S - T \ 



Our problem is to discuss these equations and to obtain 
from them expressions for the Moon's position at any time. 
The integration is best effected by observing what kinds of 
terms will disappear on substitution in the equations, and then 
assuming for the desired expressions for the coordinates a series 
of such terms multiplied by undetermined coefficients. Our 
procedure will be to discuss one by one the irregularities which 
can be isolated from one another. This will permit a survey of 
the entire field without involving needless complexity ; but if 
the Lunar Theory is to be accurate, the combinations of such 
terms with one another must also be included, and the number 
of terms employed and the labour of manipulating them becomes 
very great. 



LECTUBE III. 

THE SUN'S COOEDINATES IN TERMS OF THE TIME. 

To obtain the Moon's coordinates in terms of the time from 
the equations found in Lecture II., we must substitute in the 
expressions for the forces the developments of the Sun's co- 
ordinates which we now proceed to give. 

Employing as coordinates r', 6', of the last lecture, we have 
seen that the Sun's motion may be regarded as purely elliptical, 
so that 

a' _ 1 + e' cos (6' - ■nr') 

$'-■&' = n't -■&' + 2e' sin (n't - -&') + -. e' 2 sin 2 (n't - &') + . .. 

in which we have written for convenience n't in place of n't + e'. 
The quantities that enter the equations are 

t'\3, 



a 

~7 I > 



® 



C0S 2(6- &), 
sin 



a 'X cos ( e - 0'), 

r" / sin 

^V 008 3(0-0')- 

r / sin v ' 



Making the substitutions we find without difficulty 
fv) 3 = 1 + 1 *' 2 + W cos (n't -™' + l 



i-Y '= 1 +| e' 2 + 3e' cos (n't- -a' + \e'* cos 2 (n't -<*/) + .... 
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'' 2 (0 - 0') = ( 1 - I A C ° S 2(6- n't) 
\ 2 / sin x 



a'\ 3 cos, 
r I sin 



+ 1^ ^{2 (5 -n'() -(«'*—')} 



1 , cos 
sin 



-p^r^^-^o+c^— 0} 



+ % e' 2 ° 0S {2 (0 - n't) - 2 (n'f - «/)} 
2 sin l v 7 v /J 



+ 



• cos (0-0') = (l + 2 e ' 2 ) COS (0-^) 

r J sin ' v ' sin v 

+ 3^{(0-n'O-(n'*-«')} 

+ e 'si°n^-^ ) + ( ^- OT '^ 

+ § e' a C0S {(5 - n't) - 2 (»'* - sr')} 
8 sin c ' v /J 

+ ^ e' 2 ° 0S {(0 - n'O + 2 (»'* - W)} 
8 sin LV v " 



+ 



"'V C ° S 3 (0 - 0') = (1 - 6e' 2 ) C0S 3 (0 - n't) 
r J sin x x sin x 

+ 5e' C ° S (3 (5 - n't) - (n't - «j0) 

_ e ' C0S |3((9-n'!!) + (n'<- OT ')} 

gin i \ \ /) 

+ ^ e' 2 C .° S [3(6- n't) - 2 (n't - */)) 
8 sin l v v ' 

+ 1 e' 2 C ° S {3 (0 - n't) + 2 (n'« - w')} 

8 sin l v ' v ' 

+ 

These quantities are to be substituted where they occur in 
the expressions for the forces found in Lecture II. 

Let us now' make a few general remarks upon the result of 
the substitution. 
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It will be observed that the disturbing forces all involve the 
coefficient m'a'~ 3 . It is very important to notice that the Sun's 
parallax is not required for the evaluation of this quantity. By 
Kepler's Third Law it is derivable from observations of the 
Sun's mean motion alone. Other terms however, namely those 
with the coefficient m'/a'Hi,, involve the Sun's parallax directly; 
and that constant may be obtained by comparing the observed 
with the theoretical values of the coefficients of those in- 
equalities, with an accuracy probably greater than that of any 
other method. 

The mean disturbing force is radial, and is equal to 

1 mW 3 

2 a* { 2 e 

or the mean effect of the Sun's disturbance is to diminish the 
Moon's gravity towards the Earth ; and to diminish it more, the 
greater is the eccentricity of the Sun's orbit. Now e' has been 
diminishing for ages; hence the Moon's gravity towards the 
Earth has been increasing, and its average time for accom- 
plishing a revolution about the Earth has been diminishing. 

This is one cause of the Secular Acceleration of the Moon's 
mean motion which Halley derived from the records of ancient 
eclipses. 

It may also be noticed that the coefficient of the chief 
periodic part of the disturbing force, which involves 1—^e' 2 , 
increases as e diminishes. 

Finally let it be observed that the term with argument 
2(6- n't) + 2 (n't - w 1 ), 
which does not involve the Sun's Mean Longitude, is absent 

from the development of (— 7 ) • 2(8 — 6'). 



LECTURE IV. 

THE VABIATION. 

The Variation is the first inequality we shall consider; this 
is the inequality which is independent of eccentricities and 
mutual inclination in the orbits of the Sun and Moon. 

Let us first take the equations in the first form in which 
they are given in Lecture II., namely with t as independent 
variable : 

dP \dt) 
1 d I . d& 



r dt { dt)~ ' 



we omit the equation of motion in latitude, and in the ex- 
pressions for P, T we suppose s = 0; moreover it is possible 
and convenient to discuss separately the terms that involve 
the Sun's parallax; let these be omitted and we have 

ldV (d6V a lm' 3m' 
cM 2drd8 3m' . ... a ,. 

and if e' = 0, r' = a', m'/a' 3 = n'\ 6' = n't + e, 

\%- ®*+5-i- + I— -*<»--»-o. 

d?6 2drdd 3 ,„ . „ ,„ , 

d* + rdtdt = ~ ^sm2(d -nt-e); 

these are the equations to discuss. 
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Assume as a first approximation 

6 = nt + e + b 2 sin {2 (nt + e) - 2 (n't + e')} 
= nt + e + b 2 sin 2i/r, say ; 

- = - [1 + a 2 cos 2-Jr], 
r a r 

and we shall suppose a 2 , &2 so small that in the first instance we 
may neglect their squares and products. 

Substitute in the equations ; then 

4>(n —ri) s a 2 cos 2i|r — {v? + 4n (n — ri) 6 2 cos 2yfr\ 

+ - 3 {1 + 3a 2 cos 2^} = - m' 2 + 1 ri 2 cos 2f 



• 4 (n — ri) 2 b-i sin 2-ty + 4 (n — n') wa 2 sin 2yfr = 



3..,,.: 



n' 2 sin 2yjr. 



Hence, equating the coefficients of similar terms, we have 

¥- = n 2 4- - ri 2 
a? + 2 ' 

which gives the relation between n the Moon's mean motion, 
and - , the mean of the reciprocal of the distance ; also 



4 (n — ri 






2 —4m(n-n')b 2 = „ri* 



3 

— 4 (n — ri) 2 b 1 + 4tn(n — ri)a 2 = — n n ' 2 



•(1). 
■•(2). 



From (2) 4« (n — ri) b 2 — iria^ 
Add to (1), and substitute for fi/a?; 



3 nri 2 



2n — n" 



3 

4 (n — ri) 2 — h 2 + ^ ri 2 

3 „ 2n-ri 



3 „2n-ri 

a 2 = ~ n 2 r , 

2 n — n 



2 n — w 



11 ' 

3?i 2 - 8«ti' + 2- ri 2 



A. L. 
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- n 3 ri* 

b *- n -n' a2 + 8(n-riy 

3 ,„ n(2n-ri) 1 3 m' 2 

= — n 2 . — ? 1 — . . . + 



n' 
Calling — = m, we have 



2 — m 



a, = ~ m- 



2 '1-m" c , 11 ' 
3 — 8m + -5- m? 

7 3 2-m 1 3 m 1 



2 -(l-m) 2 ' 11 2 8(l-m) 2 ' 

v ' 3 — 8m + -=- m 2 x ' 

or, calling . , or = m,, we have 

6 n-n" 1-m 

3 , 2 + m, 

a 2 = ^m 



2'" 1 • 1 

3 — 2m! + o mi 2 



^fm^+^ + ^ + jjm, 2 . 
3 — 2m, + ~ nh 2 

These are convenient expressions, and, as it happens, very 

approximate. -If we wish to develope in ascending powers of 

m or m 1; it appears that the latter development will be the 

more convergent. 

n' 
We find by observation — = "07480, very nearly. 

Hence a 2 = -00717,95, 

o 2 = -01021,2 = 2106"-4. 

Hence the ratio of the greatest and least distances will be 

1-00717,95 : 0-99282,05, 

and the greatest angular deviation from the mean longitude 
will be 

35'6"-4, 

a very close approximation to the truth. 
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Also we have found 

!t — ™2 _L m'2 — rfl 1 1 I . 



= v? + g n' 2 = n* f 1 + Hm'j 

= m 2 x 1-00280, 
which is the relation between the actual mean motion and the 
actual mean distance (or rather mean reciprocal distance) of the 
Moon. 

Without the Sun's disturbing action, the relation between 
the mean distance and the mean motion, or rather between the 
radius of the orbit supposed circular and the uniform rate of 
angular motion along it would be 

-, = n\ 
a 6 

Hence in the actual orbit, the mean motion for a given 
mean distance is smaller than it would be without disturbance ; 

Or, for a given mean motion, the mean distance is smaller 
than it would be without disturbance. 

In fact, the relation between the mean distance and the 
mean motion is the same as it would be if the sum of the 
masses of the Earth and Moon were diminished in the ratio of 
1-00280 to 1. 



2—2 



LECTURE V. 



the variation (continued). 



We will now proceed to substitute in the differential 
equations the values of \jr and 6 which we have obtained, 
retaining terms of the order of the squares and products of 
a 2 , h 2 and m 2 or mf. 

The values to be thus substituted are 

- = - (1 + a 2 cos 2ilr), 
r a 



where 



= nt + e + 6 2 sin 2yjr, 
yfr = nt + e- {n't + e'), 
2 + rti! 



«2=2™1 



3 - 2m r + -xmr? 



Hence 



r = a 



b 2 = (l +m 1 )«, + gm 1 4 



1 — a 2 cos 2l|r + H a i (1 + cos ty) 



d 2 r 



3^ = 4a (ra — ra') a [a 2 cos 2i|r — 2a 2 2 cos 4^r], 
1 dV ri 3 

r dtf = 4 ^ ~~ n '^ 2 "** + a2 C0S 2l ^ _ 2 aa " C0S 4l ^ 
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again, 



d£ 
dt 

dd 



= n + 2 (n — n') b 2 cos 2i/r, 



-T-) = n 2 + 4tn(n- n') b 2 cos 2^r+2 (n - n'f bi [1 + cos 4i^], 



1 _ 1 

r* a? 


1 + Sa 


3 
j cos 2yfr + - ai (1 + cos 4-</r) 


Also, 


ldr 

-T7=2(n-n) 
r dt 


a 2 sin 2i|r - - a 2 2 sin 4i/r 


j 


1 dr dd ., ., 
rdidt = 2{n - n) 


na 2 sin 2ty + ■](»& — n') a 3 6 2 — 


d"0 _ 


4<(n — i 


ij b 2 sin 2i|r. 





And cos 2(0 — n't — e) = cos 2t|r — b 2 (1 — cos 4i/r), 

sin 2(0 — n't — e) = sin 2i/r + 6 2 sin 4i/r. 

Substitute these in the differential equations, and we get, 
on transposing all the terms to the left-hand sides from the first 
equation 



4 (n - n'f L a 2 2 + a 2 

-[ri i +2(n-nyb?+in(n-n')b i cos2^+2(n-nybivosty'\ 



i cos 2yfr — ^ a 2 2 cos 4t|t 



P 



3 3 

1 + „ a 2 2 + 3d2 cos 2yjr + n a£ cos 4i/r 



1 3 

- ^ rc' 2 - s m' 2 [- & 2 + COS 2l/r + 6 2 COS 4^], 



and from the second equation 
— 4 (n — w') a 6 2 s i n 2i|r 



+ 4 (m - n 1 ) 



na 2 sin 



2i|r + -j (n — n') a 2 \ — ■= naA sin 4i/r 



-I- - n' 2 [sin 2^ + 6 2 sin 4^]. 
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The coefficient of cos 2i|r in the first of these expressions, 
and that of sin 2yjr in the second, are respectively 

o 

4 (n - n'Y a 2 —4<n(n- n') 6 2 + 3 ^ a 2 - ^ w' 2 , 

3 

and — 4 (n - )i') 2 b. 2 + in (n — n') <z 2 + 5 n' 2 , 

and these are evidently reduced to zero by giving a 2 , 6 2 the 
values previously found, if we substitute for (i/aP the approxi- 
mate value m 2 + ^ m' 2 . To find the more correct value of /t/a 3 , 
equate to zero the constant term in the first expression ; 

2{n-n'fai-v?-2(n-njbi+^ (l + | aA - i re' 2 + 1 m' 2 & 2 = 0, 
that is 

g (l + |a 2 2 ) = ^ + l»'*-2(n- ? i0 2 a 2 2 + 2( w -m') 2 & 2 2 -|n' 2 6 2 



1 

2'' 



= n 2 + gw' 2 + 2(ji-m') a 



(2m 1 + m 1 2 )a a 2 -^m 1 4 



Hence we see that [i/a? differs from n 2 + ~ w' 2 only in terms 

of the fourth order, if we consider m^ a quantity of the first 
order and consequently a 2 , b 2 quantities of the second order. 
Hence also by taking 

S = W2 + ! W ' 2 ' 

in the multiplier of a 2 , when we equate to zero the coefficient 
of cos 2yjr, we only neglected a quantity of the sixth order in rn u 
and the error in the resulting values of a 2 , b 2 is of that order. 

We see that the substitution just made in our equations 
leaves outstanding terms of the fourth order in cos 4^ and 
sin 4i|r. In order to get rid of these we must add terms of this 
form to the assumed values of 1/r and 0, respectively. Suppose 
that 

- = - [1 + a 2 cos 2\fr + a 4 cos 4i|r], 
6 = nt + e + b 2 sin 2\}r + 6 4 sin 4i|r, 
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where, as we shall find, a 4 and 6 4 are small quantities of the 
fourth order. 

It may be readily seen that the additional terms introduced 
are the following : — 

1 d?r 
in - -j- 16 (n — n') 2 a 4 cos 4i|r, 

— (-n) — 8n (n — n') 6 4 cos 4i^, 

^ ^ 3a 4 cos 4<yjr, 

r 3 a? T 

d?0 

-j— - 16 (n — n'f 6 4 sin 4i/r, 

2 dr dO . ,. ... 

rdtdl 8n(»-»)a«sin4*, 

and also that we may neglect the terms added in the ex- 
pressions for 

| ra' 2 cos 2 (0 - n't - e), | n' 2 sin 2(0- »'* - e). 

If we write the terms thus produced along with the several 
terms left outstanding in our equations and then equate the 
whole to zero, we have 

1 6 (n - n'Y a 4 - 8ra (n - n') b 4 + ^ 3a 4 - 6 (ji - n'f a 2 2 - 2 (n - w') 2 &2 a 

-16(n-re') ! 64+8n(m-n')a''4+4(«-» / ) aa a 6 2- 2n ( n - w ')aa 2 

+ |w' 2 6 2 = 0, 

from which we must determine a 4 and 6 4 . 

Put ^ , = n* + I n\ 

a 3 2 

g 

i 2 = (1 + m,) a 2 + g mA 



24 LECTURES ON THE LUNAR THEORY. [LECT. V. 

and divide both equations by (n — wPf ; we get 



a t - 8 (1 + m,) bi - 6 + 2(1+ m,) 2 

27 
ai- 3 (1 +mi)m 1 2 a 2 -^mi 4 =0, 



ri6 + 3(l + 2m 1 + |m 1 2 ) 

_|(l + 2m 1 +|m 1 =) 

8(1 +m,) a 4 - 166 4 + [4 (1 + m,) - 2 (1 + m,)] a 2 2 

/ 1 \ 9 
+ 3m, 2 (l +2m 1 Ja 2 + jgm 1 4 =0. 

Simplify and multiply the last equation by ~ (1 + Mil), 

(l9 + 6m x + g ^i 2 ) a 4 - 8 (1 + m,) 6 4 - (^ + m, - ^ mA a," 

27 
- 3 (1 + mj) mj 2 a 2 - ^ m^ = 0, 

(4 + 8to! + 4m! 2 ) a 4 - 8 (1 + m,) b t + (1 + 2nh + m?) a 2 2 

+ 2 (1 + m,) ^1 + g WiJ m^ 2 a 2 + ^ (1 + ^i) ^ = 0. 

Subtract the latter from the former and 6 4 will be elimi- 
nated ; we get 

(l5 - 2m, + ^ mA a t - (— + 3m, + ^ mA a* 

- g C 1 + m ( 3 + 2 mi ) mi2cl2 ~ 32 ^ 4 + "^ mi> = °' 
which gives a t ; and this being known b 4 is found from 

&4 = 2( 1+m i) a 4+ g(l+Wh)Os 2 +Jgfl +2 m i) m i 2a 2 + 25g m i 4 - 

Taking m = '07480 as in Lecture IV, we find 
a t = -00004,580, 
& 4 = -00004,237 = 8"-740. 



LECTUKE VI. 

the variation (continued). 

Let us consider the problem of the Variation over again, 
taking now as independent variable. 

The equations of motion are given in Lecture II : — 

dhi _ jP_ _T_ du 
dd* + U ~HW Hhi 3 d0' 



„dH T 
ti 



d6 



V! 



3 ' 



, P 1r' 2 3»' ! 

where _ =M ______ cos 2(0- 0), 

so that the second equation may be written 



1 d(H>) 9 ,/dt 



= - 3n 



t\ 2 

A sin 2 (6> -ff). 



h> do \de. 

Our aim is to express t and u in terms of 6 and constant 
quantities. Now since the orbit of the Moon does not differ 
widely from a circle we may write the difference of nt + e from 
0, and the difference of au from unity as series of small periodic 
terms depending upon 6. Inspecting the form of the equations, 
it is evident that these periodic terms are of argument 2 (0 — 0') 
and its multiples ; that is 

nt + e = + periodic terms of argument 2 (0 — 0'), &c. ; 

but n't + e' = 0'; 
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therefore 

+ periodic terms of argument 2 (1 — m) — 20, &c, 
where we have written 

= e — me ; 

this constant /S is associated with (1 — m) wherever the latter 
occurs ; for brevity in writing, we shall omit it. 

We may then assume as a first approximation 

au = 1 + a 2 cos (2 — 2m) 0, 

nt + e = 0+b 2 sin(2- 2m) ; 
whence 

2(6- 6') = (2-2m)0- 2m& 2 sin (2 - 2m) 0, 

cos 2(0— 6') = mb 2 + cos (2 — 2m) 6 - m& 2 cos (4 - 4m) 0, 

sin 2 (0 - 0') = sin (2 - 2m) 6 - m& 2 sin (4 - 4m) 0, 

» * = 1 + (2 - 2m) b 2 cos (2 - 2m) 6. 

Substitute in the right-hand member of the second equa- 
tion : — 

5i ^jP- = ~ Sm2 [ sin ( 2 - 2m ) + (2 - 3m) 6 2 sin (4 - 4m) $]. 
tt do 

Therefore 

, (H*\ 3 m 2 . a . . , , 3 2 - 3m , . . - 

l0ge ("¥ j = 2 1 ^m C ° S (2 ~ 2m) 6 + 4 T^T m % C ° S (4 " 4m) °' 

which we may write 

log,, (~) = 2/t 2 cos (2 - 2m) + 2h t cos (4 - 4m) 0, 



where h is an arbitrary constant of integration, A 2 is a known 
quantity, and A 4 involves 6 2 . If we take as a second approxi- 
mation 

au = 1 + a 2 cos (2 - 2m) + a i cos (4 - 4m) 6, 

nt + e = + b 2 sin (2 - 2m) + b t sin (4 - 4m) 6, 
the above value of log* (H*/h 2 ) will not require modification 
and will supply equations of condition for determining the 
coefficients <z 2 , 6 2 , a 4 , 6 4 . 
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Thus 
so that 



dt _ n na? h 1 



log * ( w Te) = loge (t) _ I log " (?) " 2 loge aM ; 

but 

log, (n jj^J = - (1 - m) 2 6 2 a + (2 - 2m) 6 2 cos (2 - 2m) 

( + [(4 - 4m) 6 4 - (1 - m) 2 6 2 2 ] cos (4 - 4m) 6, 

a 2 / <x 2 \ 

log e aw = — -j- + a 2 cos (2 — 2m) + ( a 4 — -j- J cos (4 — 4m) 6. 

Hence we find 

-(l-m)*& 2 2 = log e (^) + i a2 2 , 

(2 — 2m) & 2 = — ^2 — 2a 2 > 
(4 - 4m) & 4 - (1 - m) 2 6 2 2 = - A 4 - 2a 4 + ^ a 2 2 . 

The remaining equations of condition that we require are 
obtained from the first equation of motion ; this may be 
written 

d- (au) 



dfr 



+ au 



Now 



1 + WwM{*n)'~*v-"\ 

d («o»( B '*y sin 2 «?-*')= g. 



d0 2 



aw = 1 4- a 2 cos (2 — 2m) + a 4 cos (4 — 4m) 0, 
whence 
a" (aw) 



d# 



• (2 - 2m) a 2 sin (2 — 2m) 6 — (4 — 4m) a 4 sin (4 — 4m) 9, 



^^ = - (2 - 2m) 2 a 2 cos (2 - 2m) 6 - (4 - 4m) 2 a 4 cos (4 - 4m) 6, 

"and 

f „' * V = m 2 [1 + (4 - 4m) b 2 cos (2 - 2m) 0], 

i' Sfsin 2(0-0')= ™ 2 [sin (2-2m)(9+(2 - 3m)Z> 2 sin(4-4m)0], 
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(n'^f cos 2 (d-0') 



= m 2 [(2 - m) b, + cos (2 - 2m) 6 + (2 - 3m) b. 2 cos (4 - 4m) 6], 
_L = I [1 + ^ _ 2A 2 cos (2 - 2m) + (A, 3 - 2A, 4 ) cos (4 - 4m) ff\. 

Substitute these in the equation above, and equate the 
coefficients of corresponding terms, 

1 + i m 2 + 1 w? (2 - m) 6 2 + | m 2 a 2 + 1 m 2 (2 - 2m) a 2 = ^(1 + /i 2 2 ) 

- (2 - 2m) 2 a 2 + (l + \ m 2 ) a 2 + 1 m 2 + (2 - 2m) m 2 6 2 = ^ (- 2A 2 ) 

- (4 - 4m) 2 a 4 + fl + ^ m 2 J a 4 + ^ m 2 (2 - 3m) 6 2 + | 'm}a 2 

- | m 2 (2 - 2m) a 2 = ^ (A 2 2 - ihj. 

If we neglect at first terms of the fourth order, we find from 
the first of these equations 

From the earlier set of equations we have 
(2 - 2m) b 2 = - h s - 2a 2 ; 
substitute this in the second equation above. We get 

- (2 - 2m) 2 + 1 +| m 2 - 2m 2 a 2 - m% + 1 m 2 = (l + g m 2 ) (- Ik,), 

or 

/ UN 3 „ n , 3 , 3 m 2 3 ,2-m 

(3-8m + y m 2 )a 2 = 1 m 2 +2^=-m 2 + 2 1 — m = 2m 2 1 ---^ > 

so that 



3 „2-m 
a 2 = h »»' 



2 l-m„ , 11 . 

u 

and 

1 1 j. 

o 2 = — , a 2 — „ — s~~ ^2 

1 — m 2 — 2m 

1 3 m 2 



1-m 2 8(l-m) 2 ' 
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These are numerically equal to the quantities denoted by 
the same symbols in Lecture IV, but 6 2 bears the contrary sign. 
We further find 

31 



(l5 - 32m + Y m *) a * 



3 rn? 

= ^ , {(1 + 3m - 2m 2 ) a 2 + (8 - 15m + 6m 2 ) &„}, 

4 1 — m ' ' '' 

or reduced to numbers 

a 4 = - -00002,210, 

6 4 = -00005,414 = 11"-17. 

Finally let us exhibit the relation between the constants 
employed in this investigation and those of Lectures IV, V; to 
distinguish them, attach accents to the latter, so that 

8 = nt + e + hi sin 2-<Jr •+- 6/ sin 4\|r, 

— = 1 +a4 cos 2i|r + a/ cos 4<y]r, 

and, omitting the constant /3 as before, 

(1 - m) = i|r + (1 - m) 6 2 ' sin 2-«|r + (1 - m) 6/ sin 4-^. 
Then 

2-f = (2 - 2m) 6 - (2 - 2m) &„' sin (2 - 2m) 6, 
sin 2-f = sin (2 - 2m) 6 - (1 - m) 6/ sin (4 - 4m) 0, 
cos 2i/r = (1 - to) b( + cos (2 - 2m) 6 - (1 - m) 6 2 ' cos (4 - 4m) 0. 

Substitute in the equation for ; we find 
nt + e = d-b 2 ' sin (2 - 2m) 6 - [6/ - (1 - m) 6 2 ' 2 ] sin (4 - 4m) 6, 
and similarly 
a'w = 1 + (1 - to) a 2 '6 2 ' + a/ cos (2 - 2m) 6 

+ [a/ - (1 — m) Oa'&jj'] cos (4 - 4m) 6. 

We observe that a' differs from a by quantities of the fourth 
order. 



LECTURE VII. 

CORRECTION OF APPROXIMATE SOLUTIONS. 

We may simplify the equations we have been dealing with, 

by a proper choice of units. Let the unit of distance be the 

radius of the circular orbit which the Moon, if undisturbed, 

would describe about the Earth in its. actual periodic time ; 

then 

/j, = nK 

Also choose the unit of time so that 

n — n' = 1, 

so that, if we take as the result of observation of the mean 

motions of the Sun and Moon, 

n':n = "07480,13, 
we get w' = "08084,9 = m 1 , 

where m, is the quantity so called in Lecture IV ; and 

fi = 1-16823,4. 
We shall frequently adopt these simplifications in what follows. 

Now let I = log e (r/a), 

so that 

rdt~dt' rdt* dt* \dt) ' r 3 a? 
and the equations discussed in Lecture IV become 



dt* + {dt)~- [dt) + a? 6 

dt 2 dt dt 

where a> = 6 — 6 . 



1 3 

2 + 2 C ° 

^ sin 2a> 



■1=0, 



= 0, 
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Now these equations are defective, for they have been 
formed by omitting certain terms from the complete equations 
as given in Lecture II. Hence, calling l , O the values of I, 0, 
which we have proved in Lecture IV to be solutions of the 
above equations, if we substitute l„, 0„ in the complete equa- 
tions of Lecture II, residuals are left, say X and Y respectively. 
And if I, 6 be solutions of the complete equations, and if we 

write 

I = k + 81, 

= d o + 80, 

where 81, S0 are small quantities whose squares and products 
may be neglected in the first instance, we obtain the following 
equations for determining SI, SO, the corrections to approximate 
solutions l , O already found : — 
„ d?Sl „ dl a dSl „ d0„ dS0 p , „,. , . 9a 

Y+ djS0 + 2 dkdS0 + 2 d0^dSl +3w 's cos 2^ = 0. 

dP dt dt dt dt 

Now let us write 

§-! + »' + .. £-£«r*-o + «, 

where c = (l + nJ+ 2 n'\ 

The quantity v consists wholly of periodic terms of the form 
cos 2iyfr multiplied by small coefficients; w contains, besides 
periodic terms, a small constant term, which however might be 
removed if we were to choose c as the constant part of fi/r s in 
place of according to the definition above. 

Let S'l, S'0 be quantities defined by the equations 
„ d*8'l _ ,. ,. dh'0 „ v , n 

„ d?8'0 ... ,.dS'l A 

then S'l, S'0 are approximations to the complete corrections SI, 
80, which if substituted in the equations that give those cor- 
rections will leave residuals, say X' and Y', where 
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x , =2 cu»cm_ 2v dve_ 3wg7 3rj/2 gin 2&)g/ 

at at oK 

T = 2 dhd^6 + 2v dM + 3 ^cos2a>8'0. 

ai at at 

We see that their value is known when h'l, h'd are determined. 
Now B'l, B'ff may be determined as follows. 
Let X = p + 2,pi cos iifr, Y=~2qismiyjr, 

where i takes all positive integral values ; and assume 
h'l = a + 2a ; cos iyjr, B'8 = 26j sin iyjr. 
Then substituting and equating coefficients, the constant 
term gives 

p - 3ca = 0, 

and the terms in iyfr give 

Pi — tfcti -2(1+ n') ibi — Sca { = 0, 
q i -i%-2(l+n')ia i = 0; 

the second of these may be written 

2(1 + r»')?-4(l + n y at -2(1 + »')i'6< = 0; 
subtract from the first and we have 

p i -2(l+ri)^-ai[i*-4 ! (l+ny + Sc] = 0, 



Pi - 2(1 +n') q f 



or a; = ■ 



3 ,: 



^-(l + nj+^rt 

and 6 i = -2(l+«')f + |- 

We see that a;, &i will be of the same order of small quantities 
as pi, qi, in general. And therefore the coefficients of the terms 
of X', Y' will be of order higher than those of X, Y. Proceed 
then to determine further corrections h"l, B"0 satisfying the 
equations 

x +- 3 r- 2 <'+»>-3r- 3cSi " - 
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then if 8'l + 8"l, 8'0 + 8"0 are substituted in the complete 
equations for 81, 80 the residuals become 

X" =2f t ^-2v^- Sw8"l + Bn'* sin 2<o8"0, 

T> = 2 M,<we ^jri 3w , 2 cog 2 

at at at 

expressions which, if developed in series of cosines and sines of 
multiples of yjr, will have coefficients of higher order than the 
corresponding coefficients in X', Y'. The like process may be 
repeated until the residuals become insensible; we then have 
sensibly correct values of 81, 80, giving 

l = l +8l = l + 8'l + 8"l + , 

= o + 80 = o + 8'0+8"0 + 

We may now take into account squares and products of the 
small quantities 81, 80 by treating l„ + 81, a + 80 as given 
approximate solutions just as we have here treated l , O ; 
substitute them in the complete equations of motion, and 
determine the residuals X, Y which they leave. These resi- 
duals will form the basis of a second approximation, and the 
operation may be repeated until no further correction is 
necessary. It is to be observed that if 81, 80 depend upon 
some such constant as the eccentricity of the Earth's orbit 
around the Sun, or the parallax of the Sun, then successive 
approximations yield correctly and separately the terms which 
depend upon the first, second, powers of that constant. 



A. L. 



LECTURE VIII. 

THE PARALLACTIC INEQUALITY. 

We shall now apply the method of the last lecture to find 
the terms in the Moon's coordinates which depend upon the 
parallax of the Sun. 

The values of I, found in Lecture IV are 
k = log,, (r/a) = - a 2 cos 2i|r , 
o = nt + e+b 2 sin 2i/r, 
and these satisfy the equations of motion in which the terms 
involving the Sun's parallax are omitted. Hence the residuals 
they leave from the complete equations are 

X = - \n'> r - || cos (0 O - 0') + M C os 3 (0„ - 6')\ , 

Y = Xn' 2 ~ || sin (0 O - 0') + ^ sin 3 (0„ - 0') 

, E-Ma 
where X = E^Ma'' 

Now from above 

O - 0' = ,}? + b 2 sin 2i/r ; 
hence we have 

sin (6 — 0') = sin yfr + = b 2 (sin yjr + sin 3^), 

COS (0 O — 6')= COS i/r — ^ 6 2 (COS i/r — COS 3i/r), 

q 

sin 3 (0„ - 0') = sin 3^ + , h (sin i/r + sin 5i/r), 

3 

COS 3 (0 O — 0') = COS 3\|r — - b 2 (cos l/r — COS Ol/r); 
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and 

^cos(0 o -0')+^cos3(0 o -0')} = (1-^6, -JJ0,,) cos* 

, /15 9 9 \ , /45 . 15 \ 

V8~ 16 2 _ 16 a 7 C0S ^ U6 2 ~ 16 ""J C0S °^' 

r - {§ sin (0„ -G') + y sin 3 (0, - 0')} = (J - ^ 6 2 - | a,) sin f 

+ (,"8 + 16 2 ~ 16 °V Sm ^ + il6 2 + 16 CT V S1U °*' 
Assume 

— SX = Xa x cos i|c + Xa 3 cos 3i/r, 
80 = Xfcj sin yjr + Xb 3 sin 3i|r, 

neglecting for the present the terms in oyjr. In the present 
case it happens that it is more advantageous to substitute 
these expressions directly in the complete equations for 81, BO 
given in the last lecture than to follow exactly the process for 
finding them by successive approximation. Omitting the 
factor \, we get 

a t cos ty + 9a 3 cos 3* + 4a 2 sin 2i|r [a x sin ijr + 3a 3 sin 3*] 

+ 3— ,[«! cos i|r+a 3 cos 3*] + 3 - s [3a 2 cos2i|r][a 1 cos T|r + a 3 cos3T^] 

-[2(1+ n') + 46 2 cos 2-f ] [6 X cos * + 36 3 cos 3i/r] 
+3n' 2 [sin 2* + b s sin 4*] [^ sin yjr + b 3 sin 3*] 

/45 15 \ 

- n' 3 I Yq b 2 — Jq as) cos 5i|r = 0, 

— &! sin i/r — 96 3 sin 3* + 4a 2 sin 2* [&! cos •* + 36 3 cos 3*] 

+ 3m' 2 [— b 2 + cos 2* + 6 2 cos 4*] [6 2 sin i|r + b a sin 3i/r] 

+[2 (1 + ri) + 46 2 cos 2i|r] [oj sin i|r + 3a 3 sin 3*] 

„ /3 21 . 3 \ . , „ /15 3 , 3 \ . „ 
+W \8 ~ ¥ 6s ~ 4 a *J sin ^ + w ( 8" + 16 h ~ T6 a V Sm 3l/r 

(45 15 \ 

— b 2 - jq cm sin 5* = 0. 

3—2 
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If we equate to zero the coefficients of cos yfr and cos 3yjt in 
the first, and those of sin yjr and sin 3i|r in the second, we obtain 
the following equations for a 1; 6 a , a 3 , 6 3 ; the terms in 5i/r remain 
outstanding, and the effect of a B , 6 5 in modifying the other co- 
efficient is neglected. 



I f 2«,-i :j^fl +|o, 



-6. 



2(l+n') + 26 k -?n'* 



6a, + 2^,0, 



+ a 3 
o I [2(l + n')-26J-6 i r 



66 2 -|«' 2 (l+6 2 )] = |^(3-96 2 -4a 2 ), 



1 + 2 n' 2 - 2a 2 + 3n% 



+a 3 [66 2 ] -bJ6a 2 ~ri*+^n'%'\ = -| w ' 2 (l-76 a -2a 2 ) ; 

s, [- 2a 2 + 1 1 a 2 J + 6, [_ 26 2 - 1 rc' 2 + | «' 2 6 

8" V" ' 2 



+ a { 9 + 3 aJ -UHl+n')} 
a, [26a] + 6, 



3 ,„/. 3, 3 \ 



2a 2 + |rc' 2 -|n' 2 6 2 



+ a,[6 (!+»')] -6 3 [9 + 3?i' 2 6 2 ] 



3... 



(s+^-^a,), 



If we require the formal values of a lt b lt a s , b s , we must 
substitute for a 2 , b 2 , p/a? the expressions we have found for 
them, and it will then be best to develope the coefficients in 
ascending powers of n'. But it is difficult to obtain by this 
process such good numerical results as we can get by sub- 
stituting the numerical values of a 2 , 6 2 , fi\a? immediately in 
the equations above. If we do so we get the equations 

4-56672a 1 -2-l72326 1 + -08093<z 3 - -051376,= §ra' s x2-87937, 

O 

2-14128« 1 -0-995646,+ -06127a 3 - -033386 s =-§jj' 2 x 0-91416, 

O 

•02349a,- •030126 1 +12-51451a 3 -6-485086 3 = |n' 2 x 500455, 

8 



02042^+ -024066,+ 6-48508a s -9000206 3 =-;Vx 5-00152. 

o 
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We notice that the first equation is not very different from 
the second doubled ; it is this fact that makes successive 
approximation a disadvantageous method and renders it ad- 
visable to include small quantities from the beginning. 

Eliminate a s , b 3 in succession from the third and fourth 
equations, thus : — ■ 

Multiply the third equation by 
9-00020h-[12-51451x 9-00020 6-48508x6-48508]= 0-127523, 
and the fourth by 

-6-48508^-[12-51451x 9-00020-648508x6 48508]=-0-091 887, 
and add ; b 3 will be eliminated. 

Again multiply the third equation by -091887 and the 
fourth by —-177317 and add; a 3 will be eliminated. Hence 
we find 

00111 9^ - -006052&! + a 3 = f n' 2 x 1-09776, • 
- -001463^ - -007034&! + b 3 = f ri* x 1-34669. 

o 

Multiply these by - "08093 and "05137 respectively and add 
to the first equation : — 

4-56672a 1 -2-l72326 1 +-08093a 3 -051376 3 =|?z' 2 x 2-87937, 

O 

-0-00009o 1 +0-000496 1 --08093a 3 =|rc' 2 x -0-08884, 

-000008a 1 -0-00036& 1 +051376 g =| n'* x 0-06918 ; 



hence 



4-56665(1! - 2-17219&! = | m' 2 x 2-85971. 

o 



Eliminate a 3} b 3 in a similar manner from the second 
equation ; 

2-14128a!-0-99o46&i+-06127a 3 --03338& 3 =-!?i' 2 x 0-91416, 

O 

-0-00007ai+0-00037&i--06127o 3 =-|n' 2 x 0-06726, 

o 

-0-00005ai-0'00023&i + -03338& 3 =-§?i' 2 x -0-04495 ; 
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hence 

214116a 1 -(V99550& 1 =-!rc' 2 x 0'93647. 
o 

From these equations we find 

a, = - 1 w' 2 x 464814 =--11392,8, 

6, = - |^x99-0.336 =--24273,4, 

a 3 = ?n' 2 x -55042= -00134,9, 

b 3 = |w' 2 x -58209= -00142,7. 



LECTURE IX. 

THE PARALLACTIC INEQUALITY (continued). 

Let us now consider the terms in 5i/r which have been left 
outstanding. 

Include additional terms \a s cos 5i/r, Xb 5 sin 5\{r in — SI, BO, 
and equate to zero the coefficients of cos 5>/r, sin 5i|r in the 
differential equations that give SI, S8. 

We have 
25a 5 + — o„ - 10 (1 + n') b 6 - 6a 2 a 3 - 6b 2 b 3 + ~ ^ aA 



-r^-^S^-iiW- ' 



- 256 6 + 10 (1+ ri) a 5 + 6a 2 6 3 + 6a 3 & 2 

+ 2 w o 3 + n ^ lfi o 2 16 
In these equations substitute 

£=l+2ra' + |«' 2 . 
Then 

(28 + 6»' + I w' 2 ) a s - 10 (1+ O 65 = (| - 9»' - ^ w' 2 ) a 2 a 3 

+ (66, + 1^)6, + ^^6,-^ a,), 
-25fc 5 flO(l+n> 5 =-6a 2 & 3 -6& 2 a 3 -| W ' 2 & 3 -«' 2 (|| h-~<>,. 
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Eliminate b* : 



a-. 



+ (7 + 2 ft >' 2 (l& 2 -~a 2 ), 
and then b 6 is given by 

From these we find 

a, = | n' a x -00595,3 = -00001 4591, 

O 

b s = ln'°- x -00710,3 ■= -00001,7410. 

o 

These numbers being so small, we see that we may safely 
ignore, as we have done, their effect in modifying the earlier 
coefficients. 

To find the effect of these coefficients upon the Moon's 

coordinates we must multiply by the factor X = p — r? ■ — . 

We shall take in accordance with the results given in 
Monthly Notices, Vol. 13, p. 177, and Appendix to the Nautical 
Almanac, 1856, 

1 = 81-5. 

Constant of Moon's Parallax = 3422"-325. 

Also we shall take in the first place, the Sun's Mean 
Parallax to be 8"'8, and in the next place 8"-9, and we will 
find the corresponding values of the coefficients of the Paral- 
lactic Inequalities. 



IX.] 
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We find 

8"-8 






8"-9 


\ = -00250,9 




\ = 


•00253,76 


\a 1 = - -00028,585 




Xa 1 = - 


-■00028,910 


\b, = - -00060,903 = - 


- 125"-62 


X6 1 = - 


- -00061,596 
= -127"05 


Xo3= -00000,3385 




Xa 3 = 


•00000,3423 


Xb s = 00000,3580 = 


0"-7384 


\6 3 = 


•00000,3620 

0"-7468 


\a s = -00000,00366 




\a 5 = 


•00000,00370 


Xb 5 = -00000,00437 = 


0" 00901 


\b s = 


•00000,00442 

0"-00911. 



These results are very fairly accurate ; but in order to get 
good values for a^, b lt we were obliged to discuss a,, b lt a 3 ,b 3 
simultaneously. Let us consider the peculiarity of the equa- 
tions from which this difficulty arose. 

Following the method of approximation of Lecture VII, if 
we neglect at first the products of 81, 80, dSl/dt, dS0/dt with the 
small quantities a 2 > b 2 , n' 2 . the equations become 

d*8l d80 y. . r 

-d¥- 2n W- 3 a° Bl + X = °> 



d*80 . dBl 

Hf +2n -dl 



+ F=0. 



Now suppose the following is a set of terms that appear 
in X pi cos (it + y), in Y q t sin (it + y), 

81 a» cos (it + 7), 88 hi sin (it + 7) ; 

then as in Lecture VII, we find 



«« = • 



Pi-Z-^i 
i 3 — n 2 + ~ n- 



7 c*n 1 

6; = -2^(1;+-^. 
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Therefore if i differs little from n, the divisor in a,{ will be 
small, and a small error or omision in the numerator of a» will 
appear magnified in the values of both a t and hi. In the case 
of the first term of the Parallactic Inequality, 

i = n — n', 



and if we take 



i 2 - n* + I n'* = - 2nn' + § n' 2 ; 



9 3 



which differ from the correct values by quantities of the fourth 
order, then 

~ n 3 ,. on - 3?/ 

and the formulae give 

3 _n'(5w-3n_0 
ai ~ 4 (n^n') (4n -!>»') ' 

, _ 2n 3 n'* 

w — n 8 (?i — w ) 2 

Now if we develope these expressions in ascending powers of m, 
i.e. n'/n, the first terms are 

15 . 15 

and these are the only terms which the formulae derived from 
our method of approximation will give correctly. 



LECTURE X. 

THE ANNUAL EQUATION. 

Let us next take into account the effect of the first power 
of the eccentricity of the Earth's orbit. We shall find that it 
produces an inequality in the Moon's coordinates, the chief part 
of which has a period of one year, and is therefore called the 
Annual Equation. 

In the formulae of Lecture VII, let the known approximate 
solutions l , d , include the Variation only ; then the equations 
for the corrections 81, S& are 

„ d?Sl . . a , dSl „ it, ._ „ . , . ,., 

lit 2 ' + 2 Sm ^ ~dt ~ s ^ 2 C ° S ^ 

- 2 [(1 + n') + 2b 2 cos 2^] ~ + 3n' 2 (sin 2i|r+ b. 2 sin 4>f)$6 = 0, 

Y+ ^ + 4a 2 sin 2f ^ + 3«' 2 (- 6, + cos 2f + & 2 cos 4>f ) 80 

+ 2[(l+n') + 2& 2 cos2^]^ =0, 

where a 2 , 6 2 , /x/a 3 are known quantities whose values are given 
in Lectures IV, V. 

Refer now to Lecture III, and we find that the terms that 
are left outstanding when the terms of the Variation are sub- 
stituted, and the parallactic terms omitted are the following : — 
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X = - 1 n'V cos (n't -n')-^- nV cos {2(0 - n't) - (n't - m/)} 

it *T 

+ | nV cos (2(0 - n't) + (n'i - w')} 

91 

F= + ^ n'V sin (2 (0 - «,'*) - (n't - u')} 

- | nV sin {2 (0 - n't) + (n'i - ot')}. 

Write a for n'i — ■&' ; then 
cos {2(6- n't) -a]= cos (2i/r - a) - (26 2 sin 2i|r) sin (2i|r - a) 
= - b 2 cos a + cos (2\jr — a) + b 2 cos (4njr — a), 
sin [2(6 — n't) — a} = + b 2 sin a -+■ sin (2i|r — a) + b 2 sin (4^ — a), 
cos (2 (6 — n't) + a] = - b 2 cos a + cos (2f + a) + b 2 cos (4i/r + a), 
sin (2 (6 - n't) + a] = - b 2 sin a + sin (2i|r + a) + b 2 sia^r + a). 

Hence / 

X = -| m ' s (l-36 2 )e'cosa-^w'Vcos(2-f-a) y \ 

+ 1 wV cos (2i/r + a)- -x- n' 2 & 2 e'cos (4i/r-a) + ^ m' 2 6/cos (4>/r +a), 

21 3 

Y = 6»i' 2 6 2 e' sin a + ^ nV sin (2-f - a) - 1 nV sin (2^ + a) 

21 3 

4- =- m' 2 & 2 e' sin (if - a) - -, n'%e' sin (if + a). 

For our present purpose we shall ignore the small terms in 
if — a and if + a which are of the sixth order. 

Assume 

— 81= a 5 e' cos a + a a e' cos (2f — a) + a 7 e' cos (2i|r + a), 

80 = 6 6 e' sin a + &„e' sin (2i|r — a) + b 7 e sin (2t|t + a). 

Now the terms which arise in the left-hand members of the 



+ 3£ 
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equations owing to terms — a p cos pt in 81, and b p sin p£ in 80, 
will be 

^) 2 ap cos pt + 2a^pa p [cos (pt — 2i|r) — cos ( pt + 2i|r)] 

-4aJ cospt + ^a 2 cos(pi-2i/r) +2a 2 cos(pi + 2i/r) 

- 2 (1 + w') p&^cos pt - 2b 2 pb p [cos (^ - 2yfr) + cos(p£ + 2i/r)] 

3 

+ s «' 2 ^ [cos (pt - 2l/f ) - cos (pt + 2l|r)], 

and 

- p% sin ptf + 2a 2 pb p [— sin (pi - 2i/r) + sin (pi + 2yjr)] 
q 

+ o W ' 2fe 2) [ Sin (P* _ 2 ^) - 2& = Sin 1 st + Sin G>* + 2 V r )] 

+ 2 (1 + n')pa p sin pi + 2b 2 pa p [sin (pi — 2i/r)+sin(pi+ 2^r)], 
respectively, neglecting the very small quantities in 4i/r. 

Hence we get the equations following : — 
Equate to zero the coefficients of cos a, sin a : 

P + ^K ~ 2n ' (! + n ') 6a + 2 (2 - n') a A + ^ . f a*. 

-2(2-n')bA + ln'%+2(2 + n')a i a 7 + m a2 € h 

- 2 (2 + n') bA + 1 n'% = | «' 2 (1 - 36 2 ) 
- v'\ - Zn'%J> t + 2 (1 + »') rc'a 6 + 2 (2 - n') a 2 6 6 

- 2 (2 - n') b 2 a e - 1 rc' 2 & 6 - 2 (2 + «') a 2 6 7 

+ 2 (2 + m') & 2 a 7 + \ n'% = - 6ri%. 
Equate the coefficients of cos (2i/r — a), sin (2i|r — a) : — 
(2 - nja s + 3 ^ a 6 - 2 (1 + n) (2 - n') b 6 + 2ra'a 2 a 5 + 3 -^ . 3 a„a, 

- (2 - n') 2 b e - Sn' 2 bA + 2(1+ »') (2 - ri) a e + 2riaA - 2n'b 2 a s 

2"'" 4 
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Equate the coefficients of cos (2i/r + a), sin (2i/r + a) : — 
(2 + nja, + S^ 3 a 7 -2(l+n')(2 + »') b 7 - 2ra'a 2 a 5 + 3 ^ . | a,a, 

(h <f," u 

- 2«'& 2 6 6 - | «,' 2 6 5 = - | « 2 , 

- (2 + nj b 7 - 3n'%b 7 + 2 (1 + n') (2 + n') a, + 2n'aA + 2w'6 2 a 5 

In equations of this class, as a general rule we would 
determine a s , b 5 approximately from the first pair, substitute 
them in the second pair and determine a„ b e approximately, 
and similarly a?, b 7 from the third pair, and repeat this approxi- 
mation as often as might be necessary. But if we refer to the 
second equation, we see that b 5 must be determined by means 
of a small divisor, and this puts any method of approximation 
at a disadvantage. In order to obtain readily satisfactory values 
for the new coefficients, we shall treat the six equations simul- 
taneously, substituting first the numerical values of the known 
quantities. 

We have found 

a., = -0071 7,95, 6 2 = 01021,20, fi,/a s = 1-17150,3. 
Hence 
3-52105a 6 - 0-l74763& 6 + -065405a, - -0293936, + -0677270? 

- -0326956 7 = | n' 2 x 0969364, 

0174763a 5 - 0-0067366 6 - -039197a, + -0177536, + -042499a, 

- -0200756 7 = | n,' 2 x - 0-040848, 
•039009a 5 + -0081536 5 + 719766a, - 4-148586, 

= |n' 2 x 3-50, 
- -001651a 5 - -0086436 6 + 414858a, -3-683356, 

= 2»i' 2 x-3-50, 
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•036687re 5 - -011455& 5 + 7-844430,-4-4981 16 7 

= |^x-0-50, 

It 

•001651a 5 + -010965&5 +4-498110,-4-330126, 

= |?i' 2 + 050. 

From the second and third pairs we find 

■016186a 6 + -007084& 5 + a e = | tj' 2 x 2-94739, 

■018678o 5 + -010326&B + 6 6 = | w' a x 426994, 
•011026o 5 - -0072036 5 + a, = | ?i' 2 x - 0-321398, 

•011072o 5 - -0100156 5 + 6 7 = | w' 2 x - 0-449338. 

Eliminate a e , b e , a 7 , b 7 from the first pair. 
Hence 

3-52015o, - -1747626 5 = |«' 2 x 0-909170, 

0-I74819o 6 - -0065366 5 =%n i x 0-003516. 
Hence 

o 5 = | «'* x - 0-70619,8 = - -00692,37, 

b 5 = | n' 2 x - 19-4268 = - -19046,3. 

Now e' is a constant found by observation; taking 
e' = 3459"-28, its value in 1850, we get 
a s e' = - -00011,61, 
6/ = _ -00319,4 = - 658"-9, 
and further 

o 6 = -03035,8, a e e'= "00050,9, 

b 6 = -04396,7, b e e'= -00073,73= 152"-09, 

c^ = - -00444,7, re/ = - -00007,457, 

b 7 = - 00623,6, & 7 e' = - 00010,46 = - 21"-57. 



LECTURE XI. 

THE EQUATION OE THE CENTRE AND THE EVECTION. 

We have seen that the equations of motion 
M 



dt* + 



in2(0-n'o] = O, 



d*0 d0dl , t 

dt* + dtdt +n 



3 

2 sir, 



are satisfied very approximately by the values 

T 

I = log - = — a 2 cos 2i|r , 

= »rf ■+ e + &2 sin 2yjr, 

where y}r = nt + e - (n't + e'), 

and a 2 , & 2 are small quantities depending upon the ratio n'jn, 
and a is a quantity depending upon n in such a way that 

while w, e are arbitrary, though subject to the assumption that 
the ratio n'/n is small. 

This solution, then, expresses a possible case of motion ; 
nevertheless it is no more than a particular case because it 
involves only two arbitrary constants, whereas the complete 
and general solution must contain four, in order that it may 
be able to satisfy any given initial conditions, that is, in order 
that the initial coordinates and their initial velocities may have 
any given values. 
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When there is no disturbance the four arbitrary constants 
are n and e, — which denote quantities similar to those expressed 
by the same symbols above, — and the two elliptic elements e 
and ot, of which e denotes the eccentricity of the orbit and 
■S7 the longitude of the apse. 

We will now shew how to complete the solution by intro- 
ducing into log (a/r) and 6 additional terms depending on 
quantities similar to e, ct, of which the former is constant and 
the latter varies slowly and uniformly with t ; and for the sake 
of simplicity we will suppose at first that e is so small that its 
square and higher powers may be neglected though it is other- 
wise arbitrary in magnitude. 

Let us assume then 

log - = a 2 cos 2i/r + e cos (nt — -as ), 

d = nt + e + b 2 sin 2i|r + 2e (1 + b ) sin (nt — ■&), 
in which the elliptic terms are of the same form as in the 
undisturbed orbit, and ■& is supposed to be slowly variable, 
so that 

dt= P > 
where p is supposed to be a small quantity of the order of the 
disturbing force. 

We will now substitute these assumed values in the differ- 
ential equations. We have 

-j- = 2(n — n) a 2 sin 2yjr + (n—p)e sin (nt — to-), 

j-„ = 4 (n — n') 2 aa cos 2^ + (n — pfe cos (nt - -a), 

dd 

-j- = n + 2 (n — ri) 6 2 cos 2i|r + 2 (n - p) (1 + b ) e cos (nt — •nr).. 

— = — 4 (n — n') b 3 sin 2i]r — 2 (n — pf (1 + b ) e sin (nt - ■&). 
ctt 

Hence 

4 (n — n'y a 2 cos 2yfr + (n — p)" e cos (nt — •sr) 

+ 2 (n — n') (n -p) <z 2 e [cos (2yfr - nt + ot) - cos (2i/r + nt — &)] 

a. l. 4 
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- {re* + 4to (n - re') b 2 cos 2^ + 4>n (to - p) (1 + b ) e cos (nt - is) 
+ 4>(n-n')(n-p)(l+b )b. 2 e[cos(2f-vt+is)+cos(2'\{r+nt— is)]} 
+ - s {1 + 3a 2 cos 2i/r} {1 + Be cos (re£ - is)} 

1 L + g cos 2f - 3 sin 2ifr [2 (1 + &o) e sin (rrf - sr)]| = 0, 



— «.'2 . 



-71- 

and 

- 4 (w - «')= 6 2 sin 2i|r - 2 («, - pf (1 + 6„) e sin (n« - «r) 
+ 4re (re - n') a 2 sin 2i|r + 2n (re - p) e sin (rei - is) 
+ 4(TO-TO')(n-p)(l+6 )ea 2 [sin(2-^-n«+OT)+sin(2i|r + ni— sr)] 
+ 2 (to - to') (re - p) e& 2 [- sin (2-v/r - re* + is)+ sin (2-f + rei - ■sr)] 

+ to' 2 L sin 2^ + 3 cos 2ifr [2 (1 + 6 ) e sin (ret - -bt)]1 = 0. 

It will of course be found that with the values of a 2 , 6 2 of 
Lecture IV, the terms independent of e vanish identically. 

Equating to zero the coefficients of cos (nt — is) in the first 
equation and sin (nt — is) in the second, we get 

(to - P y - 4to (re -p) (1 + b ) + ^ = 0, 

-2(re-p) 2 (l + & ) + 2re(re-p) =0. 

Therefore (?) — p) (1 + 6 ) = re, 

O 

and (re — pf = 4to 2 — ^ 

3 

= to 2 — ~ to' 2 , approximately, 

or - = -v m 2 = 6 > approximately. 

Now terms have been left outstanding with the arguments 
2\fr — nt-\- is, 2^ + nt — is. These may be removed by assuming 

log - = 02 cos 2ty + e cos (nt — «r) + <z 21 e cos (2i|r — wf 4- is) 

+ a^e cos (2i|r + nt — is), 

6 = nt + e + 6 2 sin 2-\fr + 2e ( 1 + b ) sin (rei — is) 

+ b. a e sin (2-v|r — ?i< + ot) + b&e sin (2i|r -f nt — is). 
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Hence in place of the former equations, we get the 
following 



(n-pY - 4m (m - p) (1 + b ) + 



Sfi 



+ 



2 (m - m') - 2w' + p) + 



2 a 3 



ClnO/tn 



2(m-m')(3m-2m'-p) + |^ 



QiaQme 



+ 



2 (n - n) (m - 2m' +p) 6 2 + 5 n' 2 



- 2 (m - m') (3w - 2m' - p) 6 2 + % m' 2 



&»=0, 



- 2 (n - pf (1 + 6.) + 2n (n - p) - 6m' 2 6 2 (1 + 6 ) 

— 2 (n — n') (m — 2n' + p) b 2 a 21 + 2(n — n') (3m — 2n' —p) b 2 a w 

2 (n - n') (n - 2m' +p) a 2 - ^ m' 2 ] b m 

+ - 2 (n -m') (3m - 2m' -p) a 2 + 1 m' 2 1 6 22 = 0. 

2m 

Multiply the second by , and add to the first ; this 

J n — p 

will eliminate 1 + b . 

, .. . „ 3/i 12mm' 2 T ., , . 

(m-p) 2 -4n 2 + -^ + 6,(1 + 6.) 

v r ' a s n—p v 7 

+ 2 (m — m') (m - 2m' +p) [a 2 a 21 — 6 2 6 21 ] 
9 u 3 

+ 2 ^ K^ + OaOaa] + jj n ' 2 &i + M 

+ 2 (m - m') (3m - 2m' - p) [a 2 a 22 - 6 2 6 22 ] 

71- 

+ 4 — — — (m - n') (m - 2n' + p) [6 2 a 21 - a 2 6 21 ] 

- 4 — — (n - ri) (3m - 2m' - p) [6 2 a 22 - a 2 6J 

+ 3m' 2 ^— [6 21 -6 22 ] = 0. 
m — p L J 

4—2 
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Also the equations obtained by equating the coefficients of 
e cos (2i|r - nt + ■&) and e sin (2yjr — nt + ■&) to zero are 

2(n-n')(n-p) + ^^a i -4 ! (n-ri)(n-p)(l+b <) )b,+3n'*(l + b ) 



+ 



(n-2n'+pf + ^ 



Oai - In (n - 2n' + p) b a = 0, 



4 (to - to') (to -p) (1 + b ) a 2 - 2 (to - to') (to -p) & 2 

- 3to' 2 (1 + b ) - (to - 2n' +pf b a + 2n(n- 2n! + p)a 21 = 0. 



Multiply the second by — 



2n 



n — 2n' + p 
this will eliminate 6 21 , and gives 



, and add to the first ; 



2(n-n')(n-p)+l^- n _ S £, + p (n-n')(n-p)(l + b )^ 



+ 



+ 



• 4 (to — to') (to - p) (1 + b ) + 4to 



to— 2?i'+p 



(n-p) 



3to' 2 + 



6to?i' 2 
to — 2to'+p 



(1 + 6.)+ 



(to-2to'+p) 2 -4to 2 + 



3/u, 



a 51 =0. 



Lastly the equations obtained by equating the coefficients of 
e cos (2i|r + nt — ra-) and e sin (2\}r + nt — ot) 
to zero are 

-2(TO-TO')(TO-|)) + |4la 2 -4(M-«')(«-p)(l + &o)&2-3TO' 2 (H-6„) 

+ (3to - 2to' -pf + 8 n a^ - 2to (3to - 2n'-p) b a = 0, 

and 

4 (to - to') (to - p) (1 + 6 ) a 2 + 2 (to - to') (to - jj) 6 2 + 3n' 2 (1 + b ) 

- (3re - 2to' -pf bn + 2to (3to - 2to' - p) a w = 0. 



Multiply the second by — , 



2to 



3to — 2n' — p 
this will eliminate 622, and gives 



and add to the first ; 
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2(„ -n')(n -J>) + JS£- 



8m 



2 a 3 3n-2n'-p 

4>n 



- 4(n - m')(ra -P)(l + &o) - 

6wn' 2 



■ 3ri* - 



3?i — 2«' — p 



3ra - 2m' - p 



(»-n')(n-p)(l+& ) 
(n — rc')(?j-p) 



+ 



(3w-2n'-p) 2 -4ii 2 + ^ 



0^ = 0. 



These six equations are to be solved by successive approxi- 
mation ; taking the first rough values of p/n and b , we find 
from the last two pairs values for a 21 , b n , a w , b^\ these are 
substituted in the first pair and yield more approximate values 
of pjn and b , and so on. 

It will be noticed that this complexity is made necessary by 
the fact that a^, b 21 are found by means of a small divisor 

(n - In' + 'pf - 4W 2 + % . 



LECTURE XII. 

THE EVECTION AND THE MOTION OF THE APSE. 

We proceed to the conversion into numbers of the formulae 
of Lecture XI. 

Take n — n' = l, 

n = 1-08084,9, 
n'= "08084,9, 

£=1-17150,3, 

a? 

log a, = 7-85609, 
log & 2 = 8-00911. 

First Approximation. 

O-p) 2 -4n 2 + 3^ = 0, 





a" 




4m 2 


4-67293,7, 




-3^ • 

a 3 


- 3-51450,9, 




(n-py = 


1-15842,8, 




n — p = 


1-07630,3, 




P = 


•00454,6, 




n — 2n' + p = 


•92369,7, 




3n-2ri-p = 


3-07630,3, 




1 + 6.= 


1 -00422,4 = «./(«■ 


-p). 
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Substitute in the equation for a 21 , 
(n - 2n' + py - 4n 2 + '-^ \a n + 2 (n - n') (11 - p) a 2 + h - 3 « 2 

Ob I Z Oj 

(n-n')(n-p)(l+b )a 2 -4<(n-'n')(n-p)(l+b )b i 



n—2n'+p 

+ 6 ""* -(l + 6b) = 0. 

n—2n +p 

The various terms give 

2(m- n')(?i— p)a 2 -01545,45 



9 u 

2^ 



•03784,83 



- 8 £-, (n - n 7 ) (» - j>) (1 + &,) a 2 - -07264,10 

- 4 (n-n')(»-p)(l +&.)&» --04415,05 

-(w-9i')(w-_p)6 3 -05144,47 



' n — 2n +p 

3n' 2 (l+6 ) -01969,25 

6 U f. (l+6 ) -04608,58 
w— 2m+jo 

•05373,43 
(w-2w' + p) 2 -85321,6 

-4n» + 3^, -1-15842,8 
a 3 

- -30521,2 



On = -17605,6. 
The equation for 6 21 is 

b * = n-2n' + p a * + 4! (n-2n'+py ( n -PK 1+b ° )a > 

- 2 (n-27 + py {n - p)k - 3 ( n-2 n n' + p^ 1+b ^ 
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2n 

o / , <hi -41201,8 

n — 2n + p 

7h ~~~ Ti 

4 (n-2n'+pf {n ~P^ 1+ 6 ») a * -03637,95 

71 — Yl! 

~ 2 r 5-7T-r.(»-p)6, --02576,42 

(n— 2n + pf r 

n' 2 

~ 3 / o /_, v (1 + 6 °) - -02308,03 

(w — 2n + pf v 7 

b a = -39955,3 

The equation for a& is 

(3n-2n' + py-4n*+3^an-2(n-n')(n-p)c h +l£a 2 

71 

~ * 3n-2n'-g {n ~ »') < » "lO &. - ^ i 1 + & ») 

-6^^ (l + 6») = 0. 

3n — 2n —p 

Here - 2 (re - n') (w - p) a 2 -'01545,45 

|^a 2 -03784,83 

- 8 o ^7 O - «') (« -p) (1 + W a. - -02181,13 

3n — 2ra— p 

_ 4 („ _ „') („ -_p) (1 + 6„) 6 2 - -04415,05 

- 4 ~ (n - n') (n -p)b 2 - "01544,69 

3n-2n -p y 

-Sn'Hl+h) --01969,25 

-6 Q ""*, -(1+ft.) -'01383,78 

- -09254,52 
(3n-2ri-p)* 9-46363,4 

-4rc 2 +3^ -115842,8 

a 3 



8-30520,6 



a„= 01114,30. 
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And 

n -n' , . . re' 2 ,.. . , . 

+ 2 (3n-2n'- P y ( n -P )b * + B ( 3n-2n'-py (1+b ° ) ' 



n 



3n-2n'-p 



•00783,011 



\sn-2n'- P )^ -^ 1+b ^ ° 0m 

8 (8»-£- f y< 1+ ^ • 00208 ' 086 

b w = -01551,37 



Second Approximation. The complete equation for n — p is 

(re -p) 2 - 4re 2 + 3 ^ 3 + 12 — & 2 (l + b ) 
v ■ r/ a 3 re-p ' 

+ 2 (re - re') (re - 2ri + p) [a 2 a 21 - W21] + k ~i [<h<hi + a 2 a^\ 
+ 4 (re - re') (re - 2re' - p) [6 2 a 21 - a 2 b 2l ] 

7b — 7) 
71 

- 4 (re - re') (3re - 2re' - p) [b 2 a^ — a 2 b^\ 

n—p 

+ 3— [&ta-6J = 0. 
re— jj 

- 4re 2 + 3 ^ - 1-15842,8 

a 3 

12 -^-6 2 (l + &„) -00080,8 

n—p 

2(re-n')(re-2re' + jp)[a 2 a 21 -&A 1 ] - -00520,1 
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9 Mr 

2 ^ L a 2«2i + OaOa] -00708,4 



3 

2 



- m' 2 [& a + 6J -00406,9 



2(«-n')(3w-2m'-^)[a 2 a 22 -& 2 6J - -00048,2 

4 ^r^("-»')(»-2n' + ii)[6 ii a ll -(i,8J - -00397,5 

_4 ^r^( n - w ')(3w-2w'-p)[6 2a22 -a 3 6 22 ] - -00003,0 

nn' 2 
3 ^r^[ 6 »~ 6 J -00756,1 

(w-p)» = 1-14859,4 
»-j) = 1-07172,5 
p= -00912,4 
p:n = -00844,2. 

Apply these numbers in the equation for 1 + 6 : — 

+ (n-p) 1 ^ ~ 2n ' ~ P) ^°" ~ "^ ~ 4 (^pf ^ 21 ~ b ^' 

1-00851,33 



n— p 



S-^^M^ + h) - -00017,51 

- ,— "L ( n - 2n ' +P) &a» - a«kJ -00086,55 
(n — pf r 

,"*""' (3n - 2»' - j>) [6,a» - a 2 & J -00000,64 
(n — pf 

-lir^p)^-^ ~ •°° 168 ' 92 

! + &„ = 1-00757,1 



Continuing the approximation for a a , b a , a.^, b w the various 
terms found are the following : — 



XII.] 



THE EVECTION AND THE MOTION OF THE APSE. 



59 



01538,88 
03784,83 
07221,53 
04410,93 
05097,33 
01975,81 
04601,14 



divisor 
■86169,5 
- 1-15842,8 
- -29673,3 



a., = -18082,0 



01538,88 
03784,83 
02182,34 
04410,93 
01540,41 
01975,82 
01390,46 



05365,53 



divisor 
9-43549,4 
- 1-15842,8 
8-27706,6 



^=•01118,03 



09254,01 



•42108,0 
•03598,79 

- -02540,22 

- -02292,94 
6 21 = -40873,6 



•00786,805 
•00328,659 
■00231,985 
•00209,403 
bn = -01556,85 



Third Approximation. We find 



1-15842,8 

- -00081,40 
■00535,02 

- -00723,54 

- -00416,02 
•00048,35 
•00407,40 
■00002,97 

- -00777,57 
(n-pf= 1-14838,0 



P 



-p = 107162,5 

p= -00922,4 

:n= -00853,5 



1-00860,75 

- -00017,57 
•00087,94 
•00000,64 

- -00167,84 
! + &„ = 1-00763,9 



01538,73 
03784,83 
07220,54 
04410,81 
05096,30 
01975,95 
04600,96 



divisor 

■86188,04 
- 1-15842,8 
- -29654,8 



05365,42 a 21 = -18092,9. 



•42128,90 
•03597,92 

- -02539,43 

- -02292,60 
6 21 = -40894,8 



01538,73 
03784,83 
02182,35 
04410,81 
01540,32 
01975,95 
01390,60 



divisor 
9-43488,0 
-1-15842,8 
8-27645,2 



- -09253,93 a a2 =-01118,10. 

•00786,881 
•00328,671 
■00231,978 
•0 0209,430 
6 22 = -01556,96 
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Fourth Approximation. 



1-15842,8 

- -00081,41 
•00534,22 

- -00727,13 

- -00416,24 
■00048,35 
•00407,66 
■00002,97 

- -00778,05 
{n-pf = 114833,2 



m-p= 1-07160,3 

p= -00924,6 

p:n= -00855,4 



1-00862,79 

- -00017,57 
■00088,04 
■00000,64 

- -0016 7,94 
1 + 6,= 1-00766,0" 



The values already found for the remaining quantities are 
sufficiently exact. 

These numbers give, taking after Hansen, 

e(l + & ) = -05491 

e = -05449 

2 (1 + 6„) e = -10982,0 = 22651"'9 

a 21 e = -00986,03 

6 21 e = -02228,44 = 4596"-6 

a 22 e = -00060,93 

b n e = -00084,85= 175"'l, 

and taking the Moon's mean annual motion 17325593", the 
annual motion of the apse is 

148202" = 41° 10' 2". 



LECTURE XIII. 

THE MOTION OP THE APSE, AND THE CHANGE OF 
THE ECCENTRICITY. 

We have seen that when the eccentricity of the Moon's 
orbit is not considered we may write 

- = -[l+a 2 cos2(0-6>')l, 
r a 

H = na 2 [l+h 2 co$2(6- &)], 

, 3 „ 2 - m 1 ,3 m 2 

where a«.= -^irv 



2 "1-m" 11 2 ' 41-m" 

3 - 8m + -g- m 2 

Let us introduce the two new arbitraries e, vr by writing 
H = hna? [1 + h 2 cos (2 - 2m) 0], 

- = jj- [1 + ciacos (2 - 2m) + e cos (6 - ot)], 

where h is a third arbitrary, which may be chosen to suit our 
convenience ; it must be unity when e = 0. 

Then 
dH dH dt 3 . /a Q ,.dt 

le = Wd0=—2 mnr - sm2 ^-^de 

= - g m 2 « 2 -g= [sin (2 - 2m) 6 + 4>me cos (2 - 2m) 8 sin (0 - w)] 

= - 1 m 2 rt 2 ^ 2 [1 - (4a 2 + A,) cos (2 - 2m) - 4e cos (6 - «■)] 
x [sin (2 - 2m) + 4me cos (2 - 2m) 6 sin (0 - «■)], 
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and also 

dd Z= d8 na2 ^ +Jh C0S ( 2 ~ 2m )^] ~na 2 h(2 - 2m) h 2 sin(2 - 2m)0. 

Now we may put h = 1 + 77, where rj vanishes with e. Neglecting 
powers of e above the first 

drj dh 



~ja = ~Th z= 3m 2 (1 + m) e sin (1 — 2md + *r) 

+ 3m 2 (1 - m) e sin (3 - 2md - w) - 9m 2 r) sin (2 - 2m) 0. 
Neglect at first the last term : 

1+wi ,--> — k~ n . ,. « „ 1— m 



77 = — 3m 2 = — jr— ecos(l — 2wi0+'st)— 3m 2 - _ ecos(3 — 2mS-w). 
1— 2m ' 3-2m • v 

Substitute this in the last term, and we get 

T3 = 3m 2 (l+m)esin(l — 2m6+'&-)+3m 2 (l— m)esin(3 — 2m6— ot) 

27 A 2 + 4m -4m 2 . .„ . 
+ T m (1- 2m) (3 -2m) 6 Sm( ^ - •> 

Now consider the other equation 

d 2 r H 2 a 1 3 

Sff = f, - $ + 2 mWr + 2 mWr C0S 2 (<? - ^ 

= -, i + n ™?n 2 r + h m 2 « 2 r [cos (2 — 2m) 6 



— 2me cos (1 — 2m0 + ra-) + 2me cos (3 — 2m# — tr)]. 

Differentiate the assumed expression for -, and let h be 

chosen so that the first differential coefficient shall have the 
same form as if h, e, ■m were constant. 

Thus 

where 

2t dh r , //-*/-» \ /n d& , n d'ur • / /\ *. « 

— T-Tfi U +a 2 cos(2 — 2m)^J + T^cos(p-OT) + e -7^ sin^t*— ?ir)=0, 
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or 
de 



d* 



-jq cos (0 — ■*?) + e tq sm(6 — ■&) = 6m 2 (1 + m) e sin (1 — 2m0 + v?) 



+ 6m 2 (1 — m)esin(3 — 2md — •&) 
2 + 4m — 4m 2 



and 



+ \ 27m i 



(1 - 2m) (3 - 2m) 



§m 3 aS e sin (0 — -as), 



dr TT H 

^ = ^(2-2m)a 2 sin(2-2m)0 + ^esin(0- OT ) 



na 



= ^ (2 - 2m) a 2 sin (2 - 2m) 



+ -r- [1 + h 2 cos (2 - 2m)0] e sin (0 - •nr), 

dV na ^ n d6 

^ = T (2-2m)^ 2 cos(2-2m)^ 

-y(2-2m)A 2 sin(2-2m)^esin(6l- OT )^ 

TO r , , ._ ,, . m , n .dd 

• + -r- [1 + A 2 cos (2- 2m) 0]e cos (0-in-W 



na, 



+ ^[l+/i 2 cos(2-2m)6>] 



_sin(0-*r)- e ^cos(0-*r) 



d0 
dt' 



Multiply by r 2 /w 2 a 3 ; then since 
,d0 



r'-=-=H = hna? [1 + h 2 cos (2 - 2m) 0], 



we have 
r 2 dV 



= (2 - 2m) 2 a 2 cos (2 - 2m) 

- (2 - 2m) A 2 sin (2 - 2m) e sin (0 - ■&) 

+ e cos (0- w) + 2A 2 cos (2 - 2m) e cos (0 — ■&) 

-™(2 - 2m) <x 2 sin(2 - 2m) 



+ [1+ A„cos(2- 2m) 0] 2 



^sin(0-w)-e^ cos(6>-ot) 
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and this is equal to 

H 2 u, 1 r" 3 r 3 

-n . » + s mS -i + a m * -, [ cos ( 2 - 2 ™) # 

nWr n 2 a 3 2 a 3 2 a 3 L v ' 



— 2me cos (1 — 2m6 + z?) + 2me cos (3 — 2m6 — ct)] 



[l+a 2 cos(2-2m)0+ecos(0+TO)] 



l + -£- + 2A 2 cos(2-2m)0 



3 



- -Tl + 5 m * C 1 - 3a * cos ( 2 - 2m)0 - 3e cos (0 - w)] 
cos (2 — 2m) + Qa 2 e cos (0 — w) 

— (■=+ 2m J ecos(l— 2m0 + -5r) — (^ — 2m J ecos(3-2m#-iB-) 

The terms in these two expressions which are independent of e 
give no new information ; equating the others : — 

- (2 - 2m) h 2 e sin (2 - 2m) sin (0 - «■) 

- r S (2 - 2m) a 2 siu (2 - 2m) 5 
h do 



+ [l + A 2 cos(2-2m)0] 2 



de . , n , <fe //, , 
sin (d—-n)-e-jn cos (0 — ot) 



d0 



3 3 

= 3m 2 ?7 — g m?e cos (0 — vr) + -? m? 



6a. 2 e cos (6 - ot) 



- (§ + 2m] e cos (1 - 2m0 + ©■) - (^ - 2mjecos(3-2m0-w) 
+ 9m 2 v cos (2 - 2m) 0. 
Reducing this expression 
_gsin(0- OT )+ e ^cos(0- OT )=(§m 2 -|m') e cos(0- OT ) 

+ (|m 2 + 3m 3 + -g-m 4 Wos(l-2m0 + 'Bj-) 

+ ( 3 m 2 - 3m 3 + -^ m 4 J e cos (3 - 2m0 - w) 
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and from before 
yq cos (0 — to) + e -jq sin (0 — to) = 18m 4 e sin (0 — to) 

+ (6m 2 + 6m 3 ) e sin (1^2m0 + to) 



+ {Qw? — 6m 3 ) e sin (3 — 2md 4 to). 
Hence, 

M= T4 W + 1 g-m^s 1 n2(^-TO) 

+ ( -^-m 2 -3m 3 - 3 rrA e sin (2 - 2m) 



+ f — -T- m 2 — g m 3 — ^ m 4 j e sin (2m# — 2to) 

^ m 2 — g m s - —~ m 4 j e sin (4 - 2m0 — 2to), 



dm 3 , 141 , / 3 „ 147 t „ , 

W = i m+ T6" m+ l 4 m - -m 4 )<:o: :^ f ,- o) 

39 



/ 9 39 \ 

+ I t m2 — ^ m3 + "o" m * ) cos (^ — % m ) & 

+ ( -j-m 2 + ^m 3 + =-^ m 4 ) cos(2m^-2TO) 



15 



+ f-^m 2 +~m 3 + ^ m 4 jcos(4-2m^-2TO). 

We notice that among these terms, one is of long period, 
approximately semiannual, and will become of greater relative 
importance than the others on integration. 

To effect this integration, assume 

II = TO + a sin 2 (0 - to)+ /3 sin (2 -2m) + y sin {(4 ~2m)0-2^}, 

so that the mean motion of II is the same as that of to, and 
substitute in the equation. 

a. l. 5 
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Then 
dU 3 . 



[lect. 



dG 



141 , 3 . 

lb 4i 



+ ^m 2 



- -=^-m 4 + 2a- Tm ! a+ -^-m 4 /3- 7 7)1=7 
x cos 2 (0 — ct) 
|m 2 - 6m 3 + ~m 4 ■)- (2 - 2m) /3 - 6m 2 a - ^ m 2 7 

x cos (2 — 2m) # 
cos(2rn0-2II) 



15 , 9 , 63 , 9 . " 



q O "I K Q Q "' 

— -^ m 2 + ^r m 3 + 3-^ m 4 + (4 — 2m) 7 — t m 2 u — ^ m 2 7 
2 2 lb 4* ^ 

x cos {(4-2m)0-2w}, 



so that if we take 



3 , 
a = — o m 2 
o 



219 4 
+ -32- m > 



_ 9 „ 15 3 99 4 
fi = ~8 m+ ~8 m 64 m ' 

3 , 3 , 51 , 

7= ^^.-^.^^ 



we have 



<ffl 3 ,309 4 , 



If we write 

this becomes 
where 



15 9 

-j- m 2 + 5 m 3 + 

m0-II = -f, 
~ = a — bcos 2-yfr, 



^m 4 lcos(2m0-2II). 



3 309 „ , 15 , 9 , , 45 . 

a = m- T m?--^m t , b = -j- m 2 + ^ m 3 + -g m 4 , 



4 "" 32 
and the solution is 



4 



tan- 1 (a /^| tan f) = Mo? ~ b * + constant. 
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Hence if we denote by -j4 the mean rate of change of ra-, we have 

QjGTt. /— , „ 

W. T7T = W - If 



dd 



3 , 225 , 4071 
:m -4 m "32 m 128' 



cfe„ 3 225 ,4071 . 

~dd= 4 m+ l2 m+ i28 m - 
We observe that a + b and a — b are the rates of separation of 
the Sun from the apse when the Sun and the apse are at 
quadratures and syzygies with one another, respectively, — that 
is if we take n for the longitude of the apse, or, what is the 
same thing, if we ignore small terms of short period. Hence 
the mean rate of separation of the Sun from the apse is a mean 
proportional between its rates when at quadratures and syzygies 
respectively with the apse*. 

[* This is the analogue for the case of the apse of Maohin and Pemberton's 
theorem on the motion of the node, inserted in the third edition of the 
Principia as a scholium to prop, xxxiii., lib. m. See some notes by Adams 
in Brewster's Life of Newton, Appendix xxx.] 
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LECTUEE XIV. 



THE LATITUDE AND THE MOTION OF THE NODE. 



Let us first treat this problem on the supposition that the 
latitude is so small that its square may be neglected. The 
equation of motion, taken from Lecture II, may be written 



d 2 z 
df'~ 



f j L + m: 



,V / 

7T{ 



1 + 



E-M r 



—, 3 cos a) 1 
r J 



E + M>, 

where z = r sin (latitude) and the cube of s is omitted; or 
neglecting the parallactic terms 



dt 



= — z 



E + VL 



The value of fi/r 3 may be considered known by the operations 
which have determined the motion in an orbit coinciding with 
the ecliptic ; that is to say, 



E = E. 
r 3 a? 



3 /3 \ 

1 + g a/ + 3a 2 cos 2i|r + ( = a£ ■+ SaA cos 4i/r 



where a has the definition of Lectures IV, V ; or numerically, 
taking 

n — n' = 1, 



£ = 117150,3 + -02523,0 cos 2i+ -00025,15 cos U. 



And 



m 



" = n' 2 = 00653,6. 



Hence 



d*z 
dP 



= -z [1-17803,9 + -02523,0 cos 2t + -00025,15 cos 4*]. 
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Let us now consider the equation 

where P = q + 2g , 1 cos It + 2g 2 cos it, 

in which q u q 2 are supposed small. 

Suppose a term in z to be c cos (kt + /3) ; when this is 
substituted in Pz there will arise terms 



c cos (k - It + /3) c cos (& + 2t + /3) 



c cos (k - 4>t + /3) c cos (A + U H- /3). 

Let us therefore assume 



z = c [cos (kt + /3) + C! cos (& + 2£ + /3) + c 2 cos (k + U + /3) + . . 



+ c_! cos (k - 2t + /3) + c_2 cos (A; - 4£ + /3) + . . .] 
c is arbitrary ; we have to determine k, c u c_ 1; &c. 

Substitute and equate coefficients: 

... + [-(k-4y+q ]c- 2 +q 1 c_ 1 +q 2 +... =0 

+q 1 C- i +[-(k-2y+q ]c- 1 + q 1 +q 2 c 1 +... = 

+q 2 C-2 +qiC-! + [-k 2 +q ]+q 1 c 1 +q 2 c 2 +... = 

+q 2 c-i + qi + [-(k+2y+q ]c 1 +q 1 c 2 +... =0 

+ q 2 + ? 1 c 1 +[-(&+4) 2 +2o]c 2 +... = 0. 

If qi, 2 , 2,--- are neglected, we have simply 

_ £• + q = ; 

this is a first approximation to the value of k. 
Taking q x into account and neglecting q 2 



c, = - 



?o- 


-(£-2) 2 ' 


2«- 


-(& + 2) 2- 
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In the actual case considered we notice that q„ does not 
differ widely from unity. Hence ft is nearly equal to unity 
also, and the denominator in c_! is small, and makes c_ x much 
more important than d. 

If we substitute these values in the third equation above, 
we have 

whence 

(ft 2 - q ) 3 - 8 (ft 2 - q o y - {16 (?. - 1) + 2 9l 2 } (ft 2 - q t ) - 8 ?1 2 = 0, 
which may be put under the form 

(ft" - q f + 2 (q - 1) (ft 2 - q t ) = - q<> + J ft « (ft 2 - J.) + g (ft 2 - ?o) s , 
whence 

^ - 1)° = (5, - iy - ?1 « - 1 ?1 2 (ft 2 -?.)+g (ft 2 - q,y. 

With this equation we can approximate very rapidly to the 
value of ft. Taking as a first approximation 

ft 2 -2» = 0, 
substitute this value of ft in the small terms and we get as a 
second approximation 

ft =1-08516,9. 
Whence the ratio of the retrograde motion of the node to the 
Moon's mean motion is 

ft/n- 1=0-1 =-00399,7, 
where g is written for ft/n. This value is very correct. Taking 
the Moon's mean annual motion as 17325593", the resulting 
annual retrograde motion of the node is 
69252" = 19° 14' 12". 
Next find the values of the coefficients c_ a , <h, c_2, c 2 . We have 

Sl = "01261,5, 
q -(k-2y= -34112,3, 
q t - (k + 2) 2 = - 8-34022,8; 
whence as a first approximation 

_, = - -03698,19, d = -00151,26. 
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Hence 

£<:_! + q a = - "00034,02, q& + & = -00014,49, 

and 

q -(k- 4) 2 = - 7-31821, q„ - (k + 4) a = - 24-6809, 

whence 

c_ 2 = - -00004,650, 

c 2 = -00000,587. 

A second approximation to c_i, c x gives 

[- (k - 2) 2 + q ] c_! = - (^ + q^ + ftc,,), 

[- (k + 2) 2 + g J d = - (ft + g 2 c_, + ftc) ; 

with the above values 

ft + q& + qiC-2 = -01261,47, ft + q$- Y + ftd = -01261,03, 

so that c_! = - -03698,00, 

e, = -00151,20. 



LECTURE XV. 

MOTION IN AN ORBIT OF ANY INCLINATION. 

Let us consider the change in the plane of the orbit 
produced in an indefinitely small time dt by the action of a 
given disturbing force. Let Z be the resolved part of the 
disturbing force at any time in a direction perpendicular to 
the plane in which the body is moving at the instant. Imagine 
the force Z to act by impulses at the small intervals of time dt, 
then Zdt will be the indefinitely small velocity generated by 
the force Z in the time dt, in the direction perpendicular to the 
plane of orbit at the instant. 

Let FP be the radius vector and P the position of the 
body at the instant. Also let PT represent the velocity at 




the instant in magnitude and direction; then if Tt be taken 
perpendicular to the plane FPT and equal to Zdt, the velocity 
and its direction after the impulse will be represented by Pt, 
and the new plane of the orbit by FPt. Draw Tm perpen- 
dicular to FP and join tm; then tmT is the angle through 
which the plane of the orbit has been turned about the radius 
vector FP in the indefinitely short time dt. 
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Now 



tmT = 



tT _ Zdt 
Tm~ v ' 

where v is the resolved part of the velocity at P perpendicular 
to the radius vector. But 

H = vr; 

hence the angle through which the orbit is turned in an 
indefinitely short time dt is 

H dt 

To find the corresponding changes in the elements that deter- 
mine the plane of the orbit, namely, the inclination of the orbit 
to a fixed plane, and the longitude of the node on that plane. 
Let NPQ be the great circle which represents the plane of the 




orbit at the time t, NR the plane of reference, usually the 
plane of the ecliptic, P the position of the body at the same 
time, and let i = PNR, the inclination, and let N be the 
longitude of the node. 

Let nPq be the position of the orbit at time t + dt. 

Take NQ = 90° ; draw nm perpendicular to NPQ and qQR 
perpendicular to NR; then QR = i; and by what we have 
just proved 



Zr 
NPn = ^. dt. 



Therefore 

where 
But 



Zr 

nm = ^fdt. sin 0, 
tL 



Zr 
qQ=-fjdt- cos 6, 



= nP. 

nm = Nn sin i = sin i dN; qQ = di. 
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Therefore 

di Zr cos 
dt = H ' 
dN _ Zr sin 
dt H sin i ' 
which give the changes of the elements required. 

Now let NMS be a spherical triangle, the centre of the 
sphere being G, the centre of gravity of the Earth and Moon ; 

M 




and let GS, GM, GN point respectively to the Sun, the Moon, 
and the node of the Moon's orbit upon the ecliptic, so that NM 
is the plane of the Moon's orbit and NS the ecliptic. Let 
M8=co, NS=0', NM= 0, of which the first is identical with 
the quantity denoted by the same symbol in Lecture II, but 
the second and third are not so. 

Then, following Lecture II, the forces on the Moon are 

£ + =£ inlf0, 

m ' r o ■ an 
^ 3 cos a>, in tsG, 

if we ignore the parallactic terms. 
This latter may be resolved into 

T - 3 cos co x (cos 8 cos ff + sin 6 sin & cos i) in MG, 

, perpendicular to 

~ 3 cos co x (sin cos & - cos sin & cos i) MG in the plane 

** of the orbit, 

, perpendicular to 

— j- 3 cos co x sin 0' sin i the plane of the 

orbit. 



XV.] MOTION IN AN ORBIT OF ANY INCLINATION. 75 

Now 

cos &) = cos cos 0' + sin sin 0' cos i 

= cos (6 - 0') cos 2 1 + cos ((9 + 00 sin 2 1 , 

sin cos 0' - cos sin 0' cos i 

= sin (0 - 0') cos 2 ~ + sin (0 + 0') sin 2 ^ . 

Hence we have the following expressions for the three forces : 

p _ /i m'r 3 m'r 
^2 "■" "J/3 2 r' 3 



T= 



£ = 



{l+cos2(0-0')}cos 4 ^ + {cos20+cos20'} 
x 2 cos 2 | sin2 | + ^ + cos 2 ^~ ^')1 sin4 | 
2 7r|_ sin 2 (0-0') cos* g 

+sin20.2cos 2 ^sin 2 |+sin2(0+0')sin 4 ^ 



3 m'r 
2 -pi 



sin t 



■sin (0-20') cos 2 ^ 

+ sin cos i + sin (0 + 20') sin 2 ^ 



Now we have seen 

di _ 2r cos 
eft - 



rfiV;__ ^rsin0 - 
eft 



H ' dt Hsini' 

also, the rate of advance of the node along the orbit is 

Zr sin 
H tan i ' 

Thus the equations of motion become 

H dd Zr sin 



together with 



r 2 eft if tan i ' 

d?r 
dl? 



d?r _m__ 



da 

d« 



= -rT. 



LECTURE XVI. 

MOTION IN AN ORBIT OF ANY INCLINATION {continued). 

To satisfy the equation at the end of Lecture XV, assume 

r = a [1+Acos 2(0-00+4,, cos 20 + A s cos W+A^os 2(0+00]. 

neglecting the square of the disturbing force and the eccen- 
tricity; thus in the small terms we write 

d6 , , dff , 

r-a, -j t =n, r=a, -^=n. 

Hence 

-~ = n?a [(2 - 2mf A l cos 2 (0 - 00 + 4>A 2 cos 20 

Cut 

+ 4m?A 3 cos 20' + (2 + 2m)M 4 cos 2 (0 + 0')] ; 
substitute in the equation 

iP-HP + r-g; 

therefore 

H* = iw, [1 + A 1 cos 2(0- 00 + A % cos 20 + A s cos 20' 

+ 4 4 cos2(0 + 0')] 
+ mVa 4 



3 

— - mWa* 



r{l + cos2(0-0')}cos 4 | 



+ {cos20+cos20'}2cos 2 |sin 2 |+{l+cos2(0 + 0')lsin 4 -2 

- n?a 4 [(2 - 2m) 2 A, cos 2(0- 0') + 4<A 2 cos 20 

+ 4wi 2 A, cos 20' + (2 + 2m) 2 A, cos 2 (0 + 0')]. 
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Again, we have the equation 

dt 
which may be written 

a 4 ["sin 2 (0 - 0') 



jrdH 3 , 
at 2 



cos* 



+ sin 20. 2 sin 2 „ cos 2 = 



+ sin 2 (0 + 0') sin 4 !]. 



Substitute for fia its approximate value w 2 a 4 in the small terms ; 
and we find from these two equations 

- (1 - m) Ai 4- 4 (1 - m) 3 A 1 + ^m 2 (l- m) cos 4 ~ = - ^ m 2 cos 4 ^ , 



■Ao + 



mi s + 



4J. 2 + 3ffl',2 cos 2 = sin 2 H 

Z Z Z 



3 i i 

= — o m 2 . 2 cos 2 - sin 2 ~ , 



3 i . i 

4m 3 J. 3 + h w^ ■ 2 cos 2 jr sin 2 „ = 0, 
2* Z Z 



- (1 + m) -4 4 + 4 (1 + m) 3 i 4 + sm ! (l + m) sin 4 1 = - ^ m 2 sin 4 1 . 

Therefore 

3 2 4 i 2-ro 

A- g m cos 2 (i_ m )(i_ 2m )(3_2r«)' 

-4 2 = — 2m 2 cos 2 5 sin 2 ^ , 

a -j 



J. 3 = 3 m 2 cos 2 ~ sin 2 



i^-rji 1 sin 4 s 



2 2(l-2m)(l + 2m)' 
2 + m 



and 

IT 2 = w 2 a 4 



2(l+wi)(l+2m)(3 + 2m)' 



3 / i *\ 

1 + m 2 - ~m 2 ( cos 4 g + siu 4 H ) 



3 m 1 



+ 5 ^ cos 4 g cos 2 (0 — 0') + 3m cos 2 5 sin 2 ~ cos 

ZiL — Til L L L 



20 



+ ■ 



3 m 2 



r ^_sin 4 -;cos2(0 + 0')J; 
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where a is defined by 



jj, — n*a 3 . 



If we preferred to define a so that the constant term in H % 
were equal to n 2 a 4 , we should have 



n*a? 



or fi = rv>a? 



= fw, + m?n 2 a 4 — = m^a" ( cos 4 -= + sin 4 g J , 



1 i 

1 + 5 m? — 3wi 2 sin 2 ~ cos 2 „ 

£i Li Li 



Let us next find the latitude and the motion of the node. 
Suppose that i = i + Ai, 

N = N tl + AN > 

in which Ai, AiV are small, i is a constant, and N varies 
slowly in proportion to the time, so that we may assume 

-df = - qn > 

AN = JVj sin 2 (5 - 0') + ^ 2 sin 20 + N s sin 20' + N* sin 2 (0 + 0'), 

Ai = A cos 2(0- 0') + 7 2 cos 20 + I s cos 20' + 1 4 cos 2 (0 + 0'). 

Then remembering that 

(W dN 

-dt =mn -dt' 

an expression that must be used in the terms of chief import- 
ance, we have 

* = ^! = _ 2 ( 1 - m) nl, sin 2 (0 - 0') - 2«J 2 sin 20 

~~ 2 ( m ~n w) nI " sin 26> ' - 2 C 1 + m ) wi " 4 sin 2 ( + 0'). 

^ = ^2 + ^ = - ^+2(l-m)ni\T 1 cos2(0-0')+2«iV 2 cos20 
etc (*C 0>C 

+ 2 (m - - ^ niV 3 cos 20'+ 2(1 + m)nN t cos 2(0+0') 
V. n at J 

in which the last term will be found to be required to get the 
constant q correctly to the order m 3 . 
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These must be equated to 

Zr cos 6 Zr sin 6 

H~' ~ Hsini 
respectively. 

Hence 

-q-2mN s i -mI a -p- 3 = --rtn?cos i; 

therefore as a first approximation 

3 

q = -jm? cos i ; 

hence 

q I 3 7T2- 2 i 

-2(l-m)/!= T ™ 2 sinicos 2 ~, 7 1 =- 5l sintcos 2 ^, 

v / 4 2 81— m 2 

3 3 

— 27, = — T m 2 sinicosi, 7 2 = -m 2 sinicosi, 

4 o 

3 m sin i 



a f 3 „ A, 3,.. r 3 mi 

2 1 m+-m 2 cos» \l s = — -rm 2 sint, J 3 = „ » 

1 + T ?«. cos i 
4 

— 2(1 + m)/ 4 = — -7m 2 sknsin 2 7r, 7 4 = s = sinisin 2 ^, 

v ' 4 2 81+m 2 

2(l-m)iV 1 = -Jm 2 cos 2 | ) ^=-|iZ^ cos2 |> 

3 3 

2i\r.j = jm 2 cosi, J\ r 2 = oncost, 

■« o 

/ 3 „ A ,, 3 „ . , T 3 racosi 



and 



2 [m+7m 2 cosi)iV3= 7m 2 cosi, iV 3 = 



8,3 

1 + -7 m cos % 

4 



2<1+*W- |m 2 sin 2 i > F i= ^sin 2 !. 

Substitute above for the quantities I 3 , N s and we get the second 
approximation to q, 

3.9 9 

q = -r ra 2 cos i — ^ m 2 cos 2 i + ttt m 3 sin 2 i. 
*• 4 32 64 
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It will be observed that J„ N 3 are of lower order than the other 
coefficients, so that in order to obtain them correctly to the 
same order as the others we were obliged to retain small terms 

in j- arising from the variability of N. 

If we take the variable plane defined by the longitude of the 
node N and the inclination i as the plane to which the position 
of the Moon is referred, we have the latitude of the Moon above 
this plane 

= Ai sin 8 — AiV sin i cos 6 



3 . . ,i 

jfflSllH COS 2 x 

o L 



1+tM cos i 
4 



1 — rn 



sin (6> -26') 



— 5 m 2 sin i cos i sin 9 

o 



+ rtll Sin I Sin 2 jr 



m 



n , 3 . 1+m 

1 + - m cos 1 
4 



sin (0 + Iff). 



LECTURE XVII. 
on hill's method of treating the lunar theory. 

Let us suppose the Moon to move in the plane of the 
ecliptic, and let us refer its motion to rectangular axes in 
rotation, the rotation being such that the axis of x passes 
always through the mean position of the Sun; that is, the 
axes rotate with angular velocity n', and if we suppose the Sun 
describes a circular orbit about the origin, its coordinates are 

x' = a, y' = 0. 

Let x, y be the coordinates of the Moon. 

Then the disturbing forces of the Sun upon the Moon 
relative to the Earth are 

m' x — a' m! m y 

parallel to the axes of x and y respectively, where 

p> = {x-aj + f, 

and the forces of the Earth on the Moon relative to the 
Earth are 

_H X _ _tu 

where r z = x l + y\ 

Now these forces may be written 

dil <m 

dx ' dy ' 

ix m' rn'x 

where i» = -H -^ . 

r p a 2 

a. l. 6 
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Bufc l-h + £ + h{"-l*) + &{*-W) + "'- 

Hence 

/* , m 7_. 1 ..^ . m 7_, 3 



We have tacitly assumed the origin to be at the centre of the 
Earth ; if we prefer to place it at the centre of gravity of the 
Earth and Moon, the necessary change is effected by multiplying 
the last terms, which correspond to the Parallactic Inequalities, 
hy {E - M)/(E + M). 

Equating these forces to the accelerations of the Moon 
parallel to the coordinate axes, we have the equations of 
motion in the form 

dt* ZU dt nX ~dx' 
d'y „ ,dx ,„ dQ, 

J + < 2n dt- n *y=dy> 

or, as they may be written, 

6?x ,dy _ dR 
d¥~ ~dt~dx~' 
d?y , dx _ dR 

dP + dt~dy~' 

where R = XI + = m' 2 (x 2 + y 2 ) 

r 2 a \ 2 " J 

Now suppose we have found values of x and y which satisfy 
this pair of equations and which involve two arbitrary constants. 
This may be accomplished by taking assumed developments 

x = 2<Zj cos i (t + 7), 

y = %bi sin *' (t + y), 
substituting in the equations, and equating coefficients of the 
various terms. The solution found will include the Variation 
and the Parallactic Inequalities. Let it be required to amend 
this solution by the introduction of the remaining two arbitrary 
constants that are required for a complete solution. 
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Let the additional terms that we seek be Bx, By, which we 
shall suppose so small that their squares and products may be 
neglected, let us consider first the terms which are multiplied 
by the first power of one of the new arbitraries, the original 
particular solution corresponding to the case in which this 
arbitrary is zero. 

Then Bx, Zy are determined by the equations 

d?Sx „ ,dSy d*R* , <PR . , „ ' 

~d^ +M dt -dxdy dX+ df dy + ' 
where X, Y are supposed known functions of x, y or of t, and 
have been added here to include disturbing causes not allowed 
for in the above form of R. 

Multiply the original equations by^7,-| and add : 

d 2 x dx d 2 y dy _dRdx dR dy 
d¥diWdi~da:~dl dy dt 
_dR 
dt ' 
since x, y are 1 the only functions of t that R involves ; whence 

'(SMS)-"* - 

where is an arbitrary constant ; this is the integral known as 
Jacobi's Integral. 

Let us write 

dx Tr . dy T7 - ■ , 
^=Fcos</>, J-Fsui*, 

then we have 7 2 = 2i? + C; 

and from the original equations themselves 



dR , dR... 

■■ -=- ,cos rf> + -j- sin <p , 
dx dy 
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. d<j) 

eft " 



(Px . 



v ~=-df sin(f> + ^ cos<l) 



dt> 



or 



(d?x dy\ . , , ld?y , _ ,dx\ 



And from these, differentiating and substituting for -=- , -^ , 
we get 



dv v dt\dt +ln J v 



d " R .A 



v d^ dVfd± \ 



+ 2 -= — =- cos <£ sin + -,-„ sin s d> 
dxdy T ay* 

-^i 81X1 ^ COS ^ 






LECTURE XVIII. 

on hill's method of treating the lunar 
theory (continued). 

The equations for Sx, Sy are 

d 2 Bx a ,dSy d*R~ d?R . 

Se-* n W = -3** B + dx!Ty^ +X ' 

tfSy . ,d$x d*R . (PR _, 
the equations for », y are 

d«« dt~ dy- 

rf /ft H'ij 
Multiply the former pair by -^-, -^ respectively, and the 

latter pair by -37-1 ~jt > an ^ add all together; we get 

dx d*8x d 2 x dhx dy d 1 Sy d 2 y dhy 
~dt ~dF + ~dP~dT + ~di dt 3 + WdJ 

_ (d?R dx d 2 R dy\ . / d>R dx d 2 R dy\ . 
~ [dx' 2 dt + dxdy dt) + [dxdy dt + dy* dt) y 

dR dSx dR d&y yd® v dy 
dx dt dy dt dt dt ' 

d?Rdx d 2 R dy _d dR 
dx' dt dxdy dt dt dx' 

d 2 R dx (PRdy _d dR 
dxdy dt dy 2 dt dt dy ' 
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Then our equation may be integrated 

dx dSx dy dSy _dR . dR . ,„ 
dt dt + dt^t~dx bx+ d^ y+1 ' 

«•-/(*£+ r$)* 

so that T is a known function of t, which involves an arbitrary 
constant. 

Now let us assume 

Sx = v cos (/> — w sin <f>, 
Sy = v sin $ + w cos <£. 

Substitute above for -=- , -f A , — =— , ,f ; we find 
a£ at dt dt 



or 



,^/cZw d&\ (dR . dR . ,\ 

v {dt- w dt) = [-d^ C0S< l )+ ^ sm V v 

I dR . , dR A _ 
+ I - -j- sm <p + -=- cos <j> J w + 1 . 

dR , dR . . dV 
But —co^^sm*^, 

dR . , dR , Tr /dd) _ , 

-^ sin ^ + ^ C0S ^ =F (i +2w . 

Therefore Fg- W f) = ^ + 7(f + 2 ,') w + , ; 



F =f- 



df eft 



„« 2«F(f + »') + T, 



whence y= \y \M +n 'l w ^ t+ \ W*^' " 

An arbitrary constant is included on the right. This 
equation shews that when w is known, v can be found ; it 
remains to determine w. 

Now by actual differentiation 



, dSx . , dSy dv dd> 

C0S *itr +sm( i > -dT = di- w -dt 



. , d'Scc , d°8y d 2 w a dv d$ (d&Y 



dt 2 ' 
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Also multiplying the differential equations for Sue, By by 
— ■ sin <j>, cos <j), respectively and adding 

. , d?8a; , d% a ,1 , dSx . . dBy\ 

- sm *iif +C0S< i > ii +2n r s *iir +sm( t , -dt) 



( d*R . , d*R A. / d*R . , d?R ,\ s 

= {-3* Ba + + to%"» V BX + {-(h : dy Sm< t , +'df C0S V ^ 



d*R 
dosdy 

— X sin $ + Y cos <f>. 



Substitute and we find 



d*w 
d? 



t dv d<f> 
'dlH, 



dt) 



^+^^1r-w\~\ +v -j^ + 2n-\-,. - w 



t dfy 
' dt* 



_,(dv 
\dt 



dt 



~S sin < f >cos< f >+ ^( cos2 ^- sin2 ^ + i sin ^ cos ^ 



d ' R ■ u 



+w 

— Xsin $ + Fcos <£. 

Now we have seen 

dv _ v dV 
di~V~dl 



dtf 

a d*R . . , , d?R " 

2 d^ Bm * 0M * + df *. 



+ 



2 {tt +n ') w + T v- 



Substitute for 2 @+ n') ^ on the left. 



We get 



dhv 

~dt' 



+ v 



dty ZdV/d± 
df + Vdt \dt 



+ w 



M^+n'l'- l'-7) Zn 4 * 



d4 
dt' 



■■■■ ■d4\ 2 
Kdt, 



dt 



+»GMr 



= v 



+w 



-2 



d*R 

ixdy 

d>R 
dxdy 



-gsin *«,.*+ £| (co^-sin^) + ^sin0cos* 

d*R . 2 , 
-d^ Sm4, 



sin <fi cos <f> 



If* 
d*R 



^^* 



-Xsin<£+ Ycos<j). 
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But by the equations proved at the end of Lecture XVII, 
the terms in v cancel one another, and we are left with the 
equation for <w : 

d?w 

-0 * " 



•©■+*S+«"-S*'* 



_ d?R . , , d*R . . 
+ 2 ^ sin * C0S ^-^ C0S * 

= -2 (Jj- + n')y.-Xsm (f>+ Ycos<j>. 

Or since 

1 dx . , 1 dy 

dd> . , 1 / dR . , dR 
dt +2n = v{-dz- Sln,}> + ^ C0S, t , 
the coefficient of w is 

3 (dRdy dR dx\ 2 _ n_ f_dRdy dR dw\ , 2 
V i \~dx~~dt + ~dydt) V* { dx dt + ~dy dt J + 



1 jdlft fdyV d?R dx dy cPR /dxV) 
~ T 1 \dar \dt) dxdy dt dt + df [dt) ] 

= P, say. 

This function P is a known function of t ; it may be seen 

that if 

x = Saj cos i (t + y), 

y = th sin i (t + 7), 

then P may be developed in the form 

P = %Ai cos i (t + 7). 

Hence if we omit the terms X, Y, due to other disturbances 
not yet allowed for, the equation for w assumes the form 

^ + w [A + A.cos (t + y)+A 2 cos 2 (t +7) + ...] = 0. 

This is identical in form with the equation treated in 
Lecture XIV, to find the motion of the node. The value of w 
may be found by the method there employed, and the value 
of v deduced from it. 
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